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1.0  INTRODUCTION  AND  SUMMARY 

Fluctuations  of  signal  and  noise  power  are  known  to 
affect  the  performance  of  passive  sonar  systems.  Prediction 
of  the  effects  of  fluctuations  are  usually  based  on  the 
distribution  of  the  signal  excess,  the  excess  of  the  signal- 
to-noise  ratio  required  at  the  receiver  input  to  achieve  a 
detection  probability  of  one-half.  If  the  signal  excess  is 
broadly  distributed,  the  detection  probability  is  approxi¬ 
mately 

co 

PD  »  J  dx  fE(x)  (1.1) 

o 

where  f£(x)  is  the  probability  density  of  the  signal  excess. 

The  formulation  above  is  based  on  the  assumption  (usually 
unstated)  that  the  relaxation  time  of  the  signal  excess 
process  is  very  large  compared  to  the  post-rectification 
averaging  time  of  the  sonar  receiver.  A  search  of  the  under¬ 
water  acoustic  literature  shows  that  in  some  cases  the  relax¬ 
ation  times  of  both  signals  and  noise  are  considerably  less 
than  a  nominal  post-rectification  averaging  time  of  five 
minutes. 

The  directive  properties  of  a  sonar  array  will  affect  the 
fluctuation  properties  of  the  noise  at  the  beamformer  output. 

If  shipping  is  the  dominant  source  of  noise,  then  the  fluct¬ 
uation  of  the  noise  envelope  at  the  beamformer  output  will  be 
greater  the  more  directive  the  array;  furthermore,  the  relax¬ 
ation  time  of  that  noise  envelope  will  be  shorter. 

This  report  develops  a  method  for  predicting  the  perfor¬ 
mance  of  a  passive  sonar  operating  with  signals  and  noise  with 
intensity  envelope  fluctuations  with  arbitrary  relaxation  times. 
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A  Monte  Carlo  method  Is  used  that  is  based  on  a  compound  random 
process  for  fluctuating  signals  and  noise,  and  on  a  multi¬ 
channel  analog. 

The  compound  random  process  for  fluctuating  noise  was 
introduced  in  Reference  1  and  is  also  described  in  Reference  2. 
This  non-Gaussian  process  appears  Gaussian  for  very  short 
observation  periods.  This  property  conforms  to  the  properties 
of  ambient  noise  as  determined  by  careful  measurement  and 
statistical  analysis  [3]. 

The  multichannel  analog,  introduced  in  Reference  1,  was 
devised  as  an  alternative  to  a  single-channel  receiving 
analog  consisting  of  an  input  bandpass  filter,  a  rectifier, 
an  averager,  and  a  fixed  threshold.  In  the  multichannel 
analog,  a  threshold  is  formed  by  a  linear  combination  of 
the  outputs  of  a  set  of  neighboring  channels.  The  analog 
is  shown  to  adapt  to  very  slow  changes  in  overall  noise 
levels  to  maintain  a  constant  false  alarm  probability.  This 
performance  characteristic  would  be  achieved  by  a  properly- 
designed  automatic  detection  system  or  by  operators  viewing 
displays  time  histories  of  multiple  inputs  such  as  beam- 
former  outputs. 

The  results  of  several  experimental  programs  concerning 
the  intensity  envelopes  of  ambient  noise  and  signals  are 
reviewed.  Depending  on  circumstances,  relaxation  times  of 
signals  or  noise  can  range  from  much  less  than  to  much  greater 
than  a  representative  post-rectification  averaging  time  of  a 
passive  sonar  receiver. 
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Methods  were  developed  for  predicting  the  probability  of 
false  alarm  for  the  case  in  which  the  inputs  have  a  common 
fluctuating  power  envelope.  An  analytical  method  is  derived 
for  predicting  the  probability  of  false  alarm  for  two  limiting 
cases:  very  slow  and  very  fast  fluctuations.  The  Monte  Carlo 
approach  is  used  for  intermediate  cases.  The  results  show 
that  the  probability  of  false  alarm  increases  with  the  rate 
of  noise  power  fluctuation  as  well  as  on  the  degree  of  fluctu¬ 
ation.  It  is  shown  that  false  alarm  probabilities  can  be  an 
order  of  magnitude  greater  than  that  predicted  assuming 
either  Gaussian  noise  or  noise  with  a  very  slowly  varying 
intensity  envelope. 

Methods  were  developed  for  predicting  the  probability 
of  detection  for  the  case  of  the  common  fluctuating  power 
envelope.  Analytical  results  are  obtained  for  the  case  of 
Gaussian  inputs,  and  for  the  case  of  very  slow  fluctuations, 
a  result  is  obtained  that  can  be  evaluated  by  numerical 
Integration.  Intermediate  cases  are  examined  by  a  Monte 
Carlo  simulation.  It  Is  shown  that  the  envelope  relaxation 
time  has  a  very  modest  effect  on  detection  differential. 
However,  the  spread  of  the  transition  curve,  which  gives  the 
probability  of  detection  as  a  function  of  input  signal-to- 
noise  power  ratio,  increases  appreciably  with  the  relaxation 
time  of  the  power  envelopes. 

The  case  in  which  the  noise  inputs  have  partially- 
correlated  noise  power  envelopes  is  examined.  Analytical 
results  are  obtained  for  the  limiting  cases  of  very  slow 
and  very  fast  fluctuations,  and  a  Monte  Carlo  simulation  is 
employed  for  intermediate  cases.  In  the  case  examined,  the 
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envelope  of  the  center  channel  had  a  correlation  of  0.75 
with  either  adjacent  channel,  one  of  which  had  a  correlation 
of  0.5  with  the  other.  For  this  case,  it  was  found  that  the 
false  alarm  probability  increased  with  the  envelope  relaxa¬ 
tion  time,  which  is  opposite  of  the  result  for  fully-correlated 
envelopes.  For  the  limiting  case  of  very  fast  envelope 
fluctuations,  the  false  alarm  probabilities  were  equal  for 
the  two  cases.  Thus,  except  for  the  limiting  case  of  very 
fast  fluctuations,  the  probability  of  false  alarm  for  the 
case  of  partially-correlated  noise  envelopes  is  greater  than 
that  for  the  case  of  fully-correlated  noise  envelopes. 

The  Monte  Carlo  approach  was  used  to  generate  transition 
curves  which  were  found  to  be  qualitatively  similar  to  those 
obtained  in  the  case  of  fully  correlated  noise  power  envelopes. 
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2.0  SIMULATION 

The  performance  of  a  passive  sonar  receiver  operating  in  an 
environment  with  highly  fluctuating  inputs  will  be  predicted 
by  simulating  the  receiver  and  its  inputs  by  means  of  digital 
computation . 

2.1  Sampled  Data  Multichannel  Analog  and  Inputs 

A  multichannel  analog  for  a  passive  sonar  receiver  was 
described  in  Ref.  1,  which  showed  that  the  analog  adapted  to 
very  slow  changes  of  noise  level  to  maintain  a  constant  false 
alarm  probability.  This  performance  characteristic  is 
achieved  by  properly-designed  automatic  detection  systems  or 
by  operators  using  displays  presenting  time  histories  of 
multiple  inputs. 

Figure  2.1  shows  one  channel  of  a  multichannel  analog 
with  two  adjacent  channels  providing  an  adaptive  offset  or 
threshold.  The  value  of  the  parameter  c  is  selected  to 
achieve  an  acceptable  false  alarm  probability.  The  number 
of  threshold  channels  need  not  be  limited  to  two.  The  symbol 
SQ  denotes  squaring,  RS  denotes  a  running  sum,  and  SI  denotes 
a  set  indicator  whose  operation  is  characterized  as  follows: 

D  «  1,  Z  >  0 

=  0,  Z.<  0  (2.1) 

The  operations  shown  in  Figure  2.2  are  equivalent  to 
those  shown  in  Figure  2.1.  The  former  are  convenient  for 
simulation  since  only  one  running  sum  is  required. 
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The  random  process  employed  for  characterizing  the 
fluctuating  inputs  was  described  in  Section  4.0  of  Ref.  1  and 
in  Ref.  2.  In  brief,  the  process  is 

N  (t )  =  VP(t)  G  (t )  (2.2) 

where  P(t)  is  a  non-negative  stationary  random  process, 

G(t)  is  zero-mean  unit  variance  stationary  Gaussian 
process  that  is  statistically  independent  of 
P(t). 

It  is  further  stipulated  that  the  relaxation  time  of  P(t), 
termed  the  power  envelope,  is  much  greater  then  that  of  G(t), 
which  then  determines  the  bandwidth  of  N(t).  The  operations 
involved  in  generating  N(t)  from  computer-generated  inputs 
are  shown  in  Figure  2.3. 

The  generation  of  a  chi-squared  envelope  is  straight 
forward;  Figure  2.4  illustrates  the  generation  of  one  with 
three  degrees  of  freedom.  The  Gaussian  inputs,  mutually 
independent,  are  applied  to  identical  linear  lowpass  filters 
(LPF).  Possibly,  the  simplest  case  is  one  in  which  these 
operations  are  running  sums  (R S).  In  that  case,  the  output 
of  the  RS  in  the  ith  channel  at  sample  time  j  is 

j 

S,(j)  =  12  G,(k)  (2.3) 

1  k=j-e+l  1 

where  G ^ ( k )  is  normal  (0,1)  and  independent  of  any  G-variate 
with  a  different  channel  subscript  and/or  time 
index 

e  is  the  number  of  terms  in  the  running  sum. 
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The  covariance  of  Sp(j)  is 


Ki(j-k) 


e  -  j+k,  0  <_  j-k  <_  e 
e  +  j-k,  e  <_  j-k  <_  0 
0,  otherwise 


(2.4) 


The  output  of  the  final  summer  is 


P(j  ) 


m 


i=l 


(2.5) 


where  m  is  the  number  of  degrees  of  freedom.  The  mean  value 
of  P(j)  is  me;  thus  e~lP(j)  is  chi-squared  with  m  degrees 
of  freedom.  The  covariance  of  P(j)  is 


Kp(j-k) 


2m  (e-j+k)2,  0  <_  j-k  _<  e 

2m  (e+j-k)2,  -e  <_  j-k  £  0 
0,  otherwise 


(2.6) 


In  the  threshold  formation  channels, 
can  be  simplified.  Figure  2.5  shows  that 
squarer  in  one  of  those  channels  is  PG2. 
squarer  to  the  G  input,  it  is  possible  to 
square-root  operation  on  the  P  input  to 
result,  as  shown  in  Figure  2.6. 


certain  operations 
the  output  of  the 
By  moving  the 
eliminate  the 
achieve  the  same 
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Further  simplifications  are  possible  for  the  case  in 
which  all  channels  have  a  common  noise  power  envelope. 

Figure  2.7  illustrates  the  input  operations  for  the  two 
threshold  formation  channels.  The  result  of  these  operations 
is  seen  to  be  a  constant  times  a  chi-squared  variate  with  two 
degrees  of  freedom.  One  half  of  the  latter  is  an  exponential 
variate  E  with  a  mean  value  of  one.  The  equivalent  output 
is  shown  on  the  right  side  of  Figure  2.8,  which  also  shows 
the  generation  of  this  output  from  the  inputs  P  and  E. 

An  exponential  variate  is  more  easily  generated  than  a 
Gaussian  variate. 

2.2  Estimation  of  Probabilities 

As  stated  previously,  the  output  D  of  the  analog  is  1 
or  0  depending  on  whether  Z  is  greater  than  zero  or  not . 

The  maximum  likelihood  estimator  of  the  probability  of  exceed¬ 
ing  zero  is 

r 

P  =  r_1  ,  (2.7) 

i=l 

and  the  indicated  sum  is  a  binomial  variate.  The  upper  95  per 
cent  confidence  limit  is  the  value  of  p  for  which 
k 

(j)pj  (l-p)r_J  =  .  025,  (2.8  ) 

J  =  0 

and  the  lower  limit  is  the  value  of  p  for  which 
r 

^(J)pJ  (l-p)r_J  =  .  025  (2.9) 


2-7 


Report  No.  4262 


Bolt  Beranek  and  Newman  Inc 


Figure  2.7  Common  Noise  Power  Envelope 


-  CPE 


Figure  2.8  Equivalent  Threshold 
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For  determing  relatively  low  values  of  false  alarm  probabilities, 
the  calculation  of  the  lower  limit  by  the  equation  above  would 
require  more  arithmetical  operations  than  necessary.  To  arrive 
at  a  more  efficient  form,  consider  the  identity 


X 


(2.10) 


J=0 

and  subtract  (2.9)  from  it: 
k-1 


£ 

j=0 


(<)pj  (l-P)r_J'  =  0.975 


(2.11) 


The  forms  (2.8)  and  (2.11)  are  more  readily  calculated  by 
recursive  methods;  for  example,  the  sum  in  (2.8)  is 

k 


.£ 


where 


s(p)  =  /-  -j  T. 

j=  0  J 

T  = 

J  j  1-p 

Tq  =  (l-p)n 


(2.12) 


J-l 


Confidence  limits  for  selected  values  of  p  are  given  in 
Table  2.1  for  r  =  1000  (independent  samples.  Values  not  given 
in  the  table  may  be  obtained  by  linear  interpolation. 


2-9 


Report  No.  4262 


Bolt  Beranek  and  Newman  Inc. 


TABLE  2.1 


Confidence  limits  for  p 


A 

P 

_Li_ 

0 

0 

.0037 

.001 

2.53-5 

.0056 

.002 

2.43-4 

.0072 

.003 

6.19-4 

.0087 

.004 

.0011 

.01 02 

.005 

.0016 

.0116 

.006 

.0022 

.0130 

.007 

.0028 

.0144 

.008 

.0035 

.0157 

.009 

.0041 

.0170 

.010 

.0048 

.0183 

.015 

.0084 

.0246 

.020 

.0123 

.0307 

.025 

.0162 

.0367 

.030 

.0203 

.0425 

.040 

.0287 

.0541 

.050 

.0373 

.0654 

.060 

.0461 

.0766 

.070 

.0550 

.0877 

.080 

.0639 

.0986 

.090 

.0730 

.1095 

.100 

.0821 

.1203 
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When  the  product  p(l-p)r  is  large  (here  p  is  the 
number  of  trials),  the  95  per  cent  confidence  limits  are 
given  approximately  by 


=  p  -  6 


=  p  +  6 


(2.13) 


(2.14) 


where 


For  p  =  0.1  and  r  =  1000,  then  p(l-p)r  =  90,  6  =  .0186, 

=  .0814,  and  pu  =  .1186.  The  corresponding  limits  from 
Table  2.1  are  p£  =  .0821  and  .1203.  Thus  for  r=  1000  and 
p  >  0.1,  the  accuracy  of  the  approximate  limits  is  satisfact- 


The  confidence  interval  pu  -  p^  is  maximum  for  p  =  0.5, 
and  it  decreases  with  a  decrease  of  p.  However,  the  ratio 
of  the  length  of  the  confidence  interval  to  the  value  of  the 
estimator  increases  with  a  decrease  of  p.  From  (2.13)  and 

(2.14)  we  obtain 


Pu"p* 


=  1.96  4)-^#- 


(2.15) 


For  r  =  1000  and  p  =  0.9,  6/p  =  .02,  and  for  f>  =  0.1,  6/p  =  19; 
thus,  for  estimating  detection  probabilities  between  0.1  and 
0.9,  sufficient  accuracy  obtains  using  1000  independent  samples. 
The  corresponding  number  of  input  samples  required  in  1000  (e+n), 
where  e  is  the  number  of  terms  in  the  running  sum  of  the 
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multichannel  analog.  If  n  =  60,  and  if  the  maximum  value  of 
e  is  60,  then  120,000  independent  input  samples  are  required 

For  the  estimation  of  false  alarm  probabilities,  the 
requirements  are  more  stringent.  For  f>  =  .01  and  r  =  1000, 
we  obtain  from  Table  2.1  6/£  =  0.67,  which  is  not  very  good, 
but  possibly  acceptable.  For  £>  =  .001,  the  lower  and  upper 
limits  are  2.5  x  10~  and  .0056  respectively,  which  for  most 
purposes  would  be  unsatisfactory.  Thus  for  the  estimation 
of  low  false  alarm  probabilities,  a  much  greater  number  of 
samples  will  be  required. 
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3.0  SPECIAL  CASE:  GAUSSIAN  NOISE  INPUTS 
3.1  Introduction 

The  special  case  of  homogeneous  stationary  Gaussian  noise 
inputs  is  of  interest  for  two  reasons:  the  results  of  an 
analysis  of  this  case  will  serve  to  check  certain  results  of 
the  simulation  and  to  show  the  convergence  of  the  output  of 
the  sampled  data  multichannel  analog  to  first-order  normality. 
The  assumption  of  normality  conditioned  on  the  power  envelope 
will  be  used  to  derive  results  in  subsequent  sections  for 
limiting  cases  of  both  slow  and  fast  fluctuations. 


Results  of  this  section  include  low-order  central 
moments,  skewness,  kurtosis,  the  probability  density  function 
of  the  output,  and  the  probability  of  false  alarm. 


For  the  special  case  of  Gaussian  noise  inputs,  the  output 
of  the  sampled-data  multichannel  analog  may  be  expressed  as 


(3.1) 


where 


is  the  average  power  of  the  input  to  channel  j , 

is  a  zero-mean  unit-variance  normal  random  variate, 
statistically  independent  of  the  other  random 
variates  with  a  different  subscript.  The  subscrict 
j  denotes  channel  numbers  and  the  subscript  i 
denotes  the  sample  frame. 


c  is  a  non-negative  constant  whose  value  affects  the 
false  alarm  probability 
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n  is  the  number  of  samples  in  the  running  sums, 

m  is  the  number  of  channels  used  for  forming  the 

threshold  for  channel  0. 


In  this  context,  the  term  homogeneous  noise  implies  that 
the  average  noise  power  in  each  channel  is  the  same.  For 
this  case,  the  output  can  be  expressed  as 

Z  =  pQ  (3-2) 

where  p  is  the  average  noise  power  of  each  channel 


Y  is  a  chi-squared  random  variable  with  n  degrees 
of  freedom, 

X  is  a  chi-squared  random  variable  with  mn  degrees 
of  freedom. 

The  probability  density  of  Z  can  be  expressed  in  terms 
of  the  probability  density  of  Q: 

fz (z)  =  P_1  fQ(zP-1 )  (3*  *0 


Furthermore,  the  moments  of  Z  can  be  expressed  in  terms  of 
the  moments  of  Q: 


miP 


(3.5) 


In  the  sequel,  the  basic  analyses  concern  the  normalized 
variable  Q,  and  the  results  are  extended  to  the  variate  Z 
via  (3.4)  or  (3*5)  as  appropriate. 
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3.2  Skewness  and  Kurtosis 

The  principal  objective  of  this  analysis  is  to  obtain  a 
quantitative  measure  of  the  non-normality  of  the  output. 

The  analysis  is  entirely  straightforward  and  is  based  on  the 
chi-squared  moments  formula  given  in  Ref.  A. 

The  mean  value  of  Q  and  the  second  through  fourth  central 
moments  are 


mQ 

n(l-c) 

(3. 

.6) 

°Q  ■ 

=  2n(14c2/'m) 

(3- 

.7) 

1 3  Q 

8n  ( 1-c 3/m2 ) 

(3- 

.8) 

'iJQ 

12n[  (l+c2/m)2n  4  i|  ( 14c1*  /m3 )  ] 

(3. 

.9) 

Note  that  the  variance  of  Q  increases  with  the  value  of  c, 
the  threshold  parameter. 

The  skewness  is  obtained  by  dividing  by  o^: 

=Q  -  2  ff  no  (3.10) 

where  f(c)  =  1  c  /'m-  3-r2 

(l+c2/m)  7 

Salient  features  of  the  function  f(c)  are  presented  in 
Table  2.1 
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TABLE  3.1 


Salient  Features  of  f(c) 


c 

-  CO 

-m 

0 

m2'3 

00 

f(c) 

)//m 

l//m+l 

1 

0 

-1A€ 

comment 

assymptotic 

. 

minimum 

maximum 

zero 

assymptotic 

I 

i 

The  case  c  =  0  corresponds  to  a  single-channel  analog;  at 
|  this  point,  the  function  f(c)  achieves  a  global  as  well  as 

a  local  maximum.  If  m  >  1,  the  function  achieves  a  global 
absolute  maximum.  Although  it  is  of  no  practical  importance, 
the  minimum  at  c  =  m  corresponds  to  an  output  with  n(m+l) 
degrees  of  freedom. 


j 

i 

i 


The  kurtosis  is  obtained  by  dividing  p^q  by  o^,  and  the 
deviation  of  the  kurtosis  from  the  value  three  for  a  normal 
random  variate  is 


6q  - 


where  g(c) 


1+c  Vm3 
( 1+c2 /m) 2 


(3.1D 


The  function  g(c)  has  even  symmetry,  and  its  salient  features 
are  shown  in  Table  2.2. 
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TABLE  3.2 

Salient  Features  of  g(c) 


c 

0 

+  m 

+  OC 

9(c) 

1 

l/(l+m) 

1/m 

comment 

maximum 

minimum 

assymptotic 

Again,  the  maximum  at  zero  is  global  for  m  >  1. 

Thus,  according  to  the  measures  of  skewness  and  kurtosis, 
the  normalized  output  of  the  multichannel  analog  is  more 
Gaussian  than  that  of  the  single-channel  analog  if  m  >  1, 
where  m  is  the  number  of  channels  employed  for  threshold 
formation.  Since  the  moments  of  the  output  Z  are  propor¬ 
tional  to  those  of  Q,  they  have  the  same  skewness  and  kurtosis. 

3.3  Probability  Density 

The  probability  density  of  the  normalized  output  will  be 
derived  from  expressions  for  its  distribution  function.  The 
distribution  of  Q  is 

FQ(q )  =  P(Q  <  q),  (3-12) 

and  substituting  (3-3)  in  3>12)  gives 

FQ(q)  =  P(Y  -  <  q)  (3-13) 
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Since  X  and  Y  are  each  derived  from  independent  sets  of 
variates,  their  joint  density  is  the  product  of  their  indi¬ 
vidual  densities. 

For  q  ^  0,  the  distribution  is 


q-r—  X 

FQ(q)  =  J  dx  fx(x)  J  dy  f y  (y )  , 


(3.14) 


and  the  probability  density  is 


fQ(q)  "  dq  FQ(q) 


w 

=  /  dx  fx(x)fy(z+|x) 


(3.15) 


If  the  appropriate  chi-squared  densities  are  substituted  in 
(3.15),  an  integral  is  obtained  tha  can  be  evaluated  by  3.331 
No.  3  of  Ref.  5  for  cases  in  which  n  is  an  even  number,  and 
the  result  obtained  is 


e^2  V  ,_c_ 

'  2  fr*  <f -:)<(£  -3-k)!  r*c 


5  -3-k 
2 


For  q  <  0,  the  distribution  is  given  by 


vq) 


dy  fy(y)  /  dx  fx(x) 


f(y-*) 


(3.17) 
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Since  X  and  Y  are  each  derived  from  independent  sets  of 
variates,  their  joint  density  is  the  product  of  their  indi¬ 
vidual  densities. 


For  q  ^  0,  the  distribution  is 

c 

00  q*« — x 

Fq  ( Q )  =  f  dx  fx(x)  J  dy  fy(y)  , 


and  the  probability  density  is 


fQ(q)  ~  Sq  FQ(q) 


co 

=  /  dx  fx(x)fy(z+£x) 


(3. 1*0 


(3-15) 


If  the  appropriate  chi-squared  densities  are  substituted  in 
(3-15),  an  integral  is  obtained  tha  can  be  evaluated  by  3*351 
No.  3  of  Ref.  5  for  cases  in  which  n  is  an  even  number,  and 
the  result  obtained  is 


Vq> 


2E  -q/2  V  [  (2  +1%  -  2-k]’ 

§  2-j  - - 

2  7T  (f  -l>u§ -wo¬ 


rn*  c 


-  -1-k 

/  <c/Dk 

t  » 

k! 


q  >  C 


For  q  £  0,  the  distribution  is  given  by 


Fs(q) 


-  00 

j  dy  fy(y)  J  dx  fx(x) 


III/  „  N 

-(y-z) 


(3-17) 
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and  the  density  is  given  by 


Vq) 


a 

I 


dy 


fy(y) 


fx[^(y-z)] 


(3-15) 


If  the  appropriate  chi-squared  densities 
(3-18),  an  integral  is  obtained  that  can 
3*351  No.  3  of  Ref.  5  for  cases  in  which 
The  result  obtained  is 


are  substituted  in 
be  evaluated  by 
n  is  an  even  number. 


-1  -k 


(3.19) 

To  the  purpose  of  calculations,  it  is  useful  to  calculate 
the  coefficients  by  means  of  recursion  formulas.  The  results 
can  be  summarized  as  follows: 


fe(q) 


(*) 


rs. 

o 


me 

2c 


£ 

k=C 


[  (l+r.)^  -2  -  k] ! 


-1 -k ) ! ( j  - 


1)! 


(«fe) 


E.a 

c  2 


Pi,(q)  .  q  <  o 


k=Ci 


Pk(0)  ,  q  »  0 
n  1 
-q/2 

e  2^  Qk(q)  >  q  >  0 


(3-20) 


where 


PkU) 


C-f  -  k) 


[ (l+m)§  -1  -  k]k 


(i&s)  ^Pk-l(q)  (3.2D 


Vq) 


(^-  k) 


-1 


[ (l+m)|  -1 -k  ]k 


&  f  «k-l(q)  (3-22) 


n-l 
r7  1 


P0(q)  -  Q0(q)  -  f 


TTm 


mn 

2 


2  1 


(ST?)  fe)  (3.23) 


The^c  results  are  valid  only  when  n  is  an  even  number. 
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Figure  3.1  through  3-^  show  the  density  of  Q  for  various 
values  of  n,  the  number  of  terms  in  the  running  sum,  and  c, 
the  threshold  parameter.  The  number  m  of  threshold  channels 
is  two.  To  show  the  departure  from  normality,  normal  curves 
with  the  same  mean  and  variance  are  also  plotted.  These 
functions  are  given  by 


fA(z) 


Z 


z+n ( c-3  ) 
2n  ( 1  +  c 2  /rr. 


) 


^ 2 n ( 1+c2 /m) 


(3.2*0 


where  Z(  )  is  the  normal  density  function  for  zero  mean  and 
unit  variance. 


Table  3.3  gives  low-order  moments  and  measures  for  the 
same  combinations  of  parameter  values.  Even  for  n  =  6,  the 
deviation  from  normal  is  not  extreme;  however,  the  convergence 
is  slow,  and  the  deviation  is  still  apparent  for  n  =  60. 


TABLE  3.3 

Measures  and  Moments 


n 

6 

6 

14 

14 

30 

30 

60 

60 

c 

1.0 

1.5 

1.0 

1.4 

1.0 

1.3 

1.0 

1.2 

mQ 

0 

-3.0 

0 

-5.6 

0 

-9.0 

0 

-12.0 

°Q 

4.24 

5.05 

6.48 

7.45 

9.49 

10.5 

13.4 

14.4 

*Q 

0.47 

.06 

0.31 

.09 

0.21 

.09 

.15 

.09 

«0 

1.0 

.72 

0.43 

.32 

0.20 

.16 

.10 

.09 
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3.4  False  Alarm  Probability 

The  false  alarm  probability  is 

pf  «  P(Z  >  ojs) 

00 

«  J  d2  f2(z|S)  (3.25) 

where  S  denotes  that  no  signal  is  present. 

Substituting  (3.*0  in  (3.25)  gives 

CD 

FF  *  p'1  /dz  fq(zp-‘ )  (3. 20 

0 

A  change  of  variable  x  =  zp-1  gives 

oc 

Fp  *  £  dzfj-Cz)  ,  (3-2‘,' 

which  shows  that  the  false  alarm  probability  is  independent 
of  the  average  input  power. 


Substituting  (3.16)  in  (3.27)  yields  a  form  that  can  be 
integrated  using  3-351  No.  2  of  Ref.  5-  Then  if  the  summation 
index  is  changed  to  j  =  n/2-l-k  the  result  is 


nn  2  (Hu  - 1  +  k )’ 
p  ,  V  I  c  \ 

*  h  ^ 

T_1+k  c 

RJ  k  m+c  RJ-1  *  J  -  1 


im 

\rr*c  / 


(3.2°) 


(3-30) 

(3.3D 
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If  the  analog  output  is  assumed  to  be  normal,  then  the 
false  alarm  probability  is 


=  Q 


n(c-l) 


l+cVm) 


where  Q(x) 


Z(u) 


(3.32) 


Figures  3.5  through  3.8  show  the  false  alarm  probability 
of  the  sampled-data  multichannel  analog  as  a  function  of  the 
threshold  parameter  C.  Also  shown  is  the  false  alarm  prob¬ 
ability  calculated  assuming  a  normal  output.  The  sequence 
shows,  as  expected,  that  the  approximation  improves  with  n,  the 
number  of  terms  in  the  past-rectification  running  sum.  Or, 
put  another  way,  it  shows  the  convergence  of  the  output  density 
to  normal. 
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4.0  NOISE  AND  SIGNAL  CHARACTERIZATION 
4.1  Introduction 

The  results  of  numerous  experiments  show  that  underwater 
ambient  (acoustic)  noise  cannot  be  represented  by  a  stationary 
Gaussian  process.  Most  of  these  experiments  concern  the  squared 
pressure  averaged  over  a  period  of  time  on  the  order  of  a 
minute.  The  excessive  variance  of  that  quantity  (or  ten  times 
the  logarithm  of  that  quantity  to  the  base  ten)  shows  that 
the  noise  pressure  does  not  have  stationary  Gaussian  character¬ 
istics.  This  conclusion  is  also  supported  by  a  much  smaller 
number  of  experiments  that  concern  the  statistics  of  the  noise 
pressure  itself.  These  experiments  provide  the  additional 
information  that  the  noise  pressure  exhibits  nearly  Gaussian 
characteristics  for  short  observation  periods. 

A  stochastic  process  displaying  the  characteristics 
described  above  is  easily  evolved.  First  define  G(t)  as  a 
stationary  Gaussian  process  with  zero  mean  and  unit  variance. 

Next  multiply  it  by  /v,  v>0;  the  result  is  a  stationary  Gaussian 
process  with  zero  mean  and  variance  v.  Instead,  multiply  G(t) 
by  /V,  where  V  is  a  non-negative  random  variable  statistically 
independent  of  G(t).  This  is  a  zero-mean  non-Gaussian  process 
called  a  compound  Gaussian  process  that  has  been  studied  and 
applied  [6,7].  Given  a  realization  of  V,  the  process  is  condi¬ 
tionally  a  stationary  zero-mean  Gaussian  process.  This  process 
does  not  have  the  characteristics  described  in  the  previous 
paragraph  because  a  sample  function  of  the  process  has  stationary 
Gaussian  characteristics  regardless  of  the  period  of  observation. 
For  the  final  step  of  the  evolution,  multiply  G(t)  by  /V ( t  ) , 
where  V(t)  is  a  non-negative  random  process  independent  of  G(t). 
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The  result  Is  again  a  zero-mean  non-Gaussian  process.  It  is 
stationary  if  V(t)  is  stationary,  and  vice  versa.  Now  suppose 
that  an  observation  period  Tq  can  be  found  such  that 


V 


>>  T„  >>  T, 


(4.1) 


where  Ty  is  the  relaxation  time  of  the  V-process 
is  the  relaxation  time  of  the  G-process. 

U 

Then  during  the  observation  period,  a  sample  function  of  the 
V-process  will  vary  hardly  at  all,  and  a  sample  function  of 
the  G-process  will  go  through  many  alternations.  Thus,  the 
sample  function  will  exhibit  nearly  Gaussian  characteristics 
during  a  period  of  observation  period  satisfying  (4.1).  Arase 
and  Arase  [3]  found  that  sample  functions  of  ambient  noise 
pressure  satisfied  stationary  Gaussian  criteria  more  often 
with  10-second  observation  periods  than  with  periods  of  20 
and  40  seconds. 


For  the  most  general  representation,  the  V-process  would 
be  nonstationary  because  the  phenomena  that  generate  and  propa¬ 
gate  sound  may  have  annual,  lunar,  or  diurnal  dependencies. 
However,  if  the  period  of  interest  is  very  small  compared  to 
the  period  of  the  shortest  cycle,  then  the  process,  with 
suitable  selection  of  parameters,  can  be  assumed  to  be  locally 
stationary.  For  example,  the  post-rectification  averaging 
time  of  a  passive  sonar  is  very  small  compared  to  this  diurnal 
period . 
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The  probability  density  function  for  the  compound  process 
N  ( t )  =  /Vlt7  G  (t  )* 

may  be  found  from  its  conditional  density  given  V(t),  which 
is  Gaussian: 

fN  [x  |  V  ( t )  =  v  ( t )  ]  =  1  exp  -  i  /  r 

N  /  27iv(t)  2  V  ( t )  '  ' 2 } 

The  joint  density  of  N(t)  and  V(t)  is 

fNV(x,v;t)  =  fN[x|V(t)  =  v(t)]  f v  ( v ;  t )  ,  (4.3) 

and  the  marginal  density  of  N(t)  is 

fN(x;t)  =  J  dv  f v ( v ; t )  fN[x|V(t)  =  v(t)]  (4.4) 

•Cod 

Substituting  (4.2)  in  (4.4)  gives 

00 

fN(x;t)  =  J  dv  fr(v.t)  exp  -l  ^  (4.5) 

—  00  ^ 

A  change  of  variable  v  =  y2  gives  the  alternative  form  (4.6)  of 

[1]. 


The  compound  process  N(t)  =  /V ( t )  G(t)  has  been  described 
previously  [1,2].  Another  stochastic  process  has  been  proposed 
[8]  for  arctic  ambient  noise:  /v  G(t)  +  I(t),  where  the  first 
term  is  stationary  Gaussian  noise,  zero  mean  and  variance  V, 
and  I(t)  is  an  impulsive  monstationary  stochastic  process. 

In  the  cited  reference,  the  statistical  properties  of  the 

*  In  other  sections  of  this  report,  V(t)  is  denoted  by  P(tl. 
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proposed  process  were  not  elaborated.  Certain  statistical 
properties  of  N(t),  previously  derived,  are  summarized  below. 


The  first-order  moments  can  be  obtained  by  substituting 
(4.5)  in 

miN(t)  =  y*dx  fN(x;t)  ,  (4.0 

or,  perhaps  more  easily  from 

m.N(t)  =  E[Vi/2(t)]  E]Gi(t)].  (4.7) 


The  moments  of  the  standardized  Gaussian  variate  are  well 
known;  thus 


miN(t) 


•  0,  i 


=  m 


odd 

i! 


|V  (|)!  2i/2 


i  even 


(4.8) 


The  mean  value  is  zero,  as  would  be  expected  for  either  a 
bandpass  or  highpass  process;  this  was  also  demonstrated 
experimentally  [9].  Since  the  mean  value  is  zero,  all  moments 
are  central  moments.  The  skewness  coefficient  s,.,  defined  as 
U3N  t  a*  ,  is  zero. 

The  kurtosis  is 


m2  +  o2 


3  (1  +  — ) 

m2 

v 


(4.9) 
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The  values  of  sM  and  k>7  are  consistent  with  experimental  results 
reported  in  (3)  and  discussed  in  (1  and  2). 


Consider  the  second-order  properties  of  N(t)  for  the  cases 
in  which  V(t)  is  a  stationary  process.  The  covariance  of  K(t) 
is 


Kn(t)  =  E  [/  V(t )  V(t+tT]  E  CG(t)  G(t+x )]  . 

(h.io) 

If  Tv 
with  x 

>>  To  ,  then  the  second  factor  decays  much  more 

(j 

than  does  the  first,  and 

rapidly 

Kj^(x)  »  m\’°G  ( 1  )  > 

(A. 11) 

where 

m  =  E[  V  ( t )  ] 

V 

Pq(t)  =  E [ G  ( t )  G  (t  +  x  )  ] 

Now  let  Q(t)  =  N  (t);  then 

Kq(t)  A  E  [Q(t )  Q(t+x)] 

*  [1  +  2p*(x)3  Ky(x)  +  2[iry?G(x)] 

(4.12) 

~  Kv(x)  ,  x  »  TG 

(4.13) 
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Note  that  the  presence  of  the  V-process  is  evident  from  K^(t), 
but  not  from  Kn(t),  given  Ty  >>  TQ. 

If  an  estimate  of  K^(-r)  is  derived  from  measurements,  then 
from  (4.11) 

fny  =  Kn(0)  (4.14) 

Pg(t)  =  Kn(t)/Kn(0)  (4.15) 

And  if  an  estimate  of  K^(  )  is  also  obtained,  then  use  of 
(4.12),  (4.14),  and  (4.15)  yields 


Kv(t) 


ij(0,  V1L-JV2 

K'(0)  ♦  28*(t> 


(H.lf  ) 


That  Is  the  pair  of  covariance  functions  K^(t)  and  K^(t) 
specifies  the  component  process  V(t)  and  G(t)  to  the  second 
order,  given  Tv  >>  Tp. 

To  this  point,  only  two  specifications  were  placed  on  the 
V-process:  V(t)  >  0,  and  Tv  >>  TQ .  A  normalized  chi-squarec 

process  has  been  proposed  for  the  V-process  [1,2]: 

n 

V(t)  =  mvn_1EGi(t)  (4.17) 

i  =  l 

where  >  0 

Gi(t)  is  stationary  Gauss  [O.r(T)],  r(0)  =  1.  independent 
of  G (t )  and  Gj  (t ) ,  i  ¥  j 
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The  covariance  function  of  V(t)  is 


Ky( t )  =  2m^n_ 1 r2 ( r )  (4.18) 

Given  that  an  estimate  of  Ky(  )  is  available  as  discussed  in  the 
previous  paragraph,  then  ( h . 1 8 )  and  (4.14)  give 


n  =  2Kj(0)/Kv(0)  (4.18) 

from  which  the  nearest  integer  value  would  be  selected.  And 
substituting  (4.19)  and  (4.14)  in  (4.18)  and  solving  gives 

r2 ( r )  =  Kv(t)/Kv(0)  (4.20) 

Equations  (4.14)  through  (4.20)  show  that  the  pair  of  factor 
processes  are  completely  specified  by  the  pair  of  covariance 
functions  K^(t)  and  K^(t).  If  N(t)  were  a  Gaussian  process, 
only  K^(t)  would  be  required  for  its  specification. 

A  search  of  the  literature  has  not  revealed  any  joint 
estimates  of  K^(t)  and  K^(t). 


4.2  Envelope  Estimation 

The  square  of  the  compound  process  is 

Q(t)  &  N2 (t)  =  V(t)  G2 (t)  (4.21) 

Given  that  Ty  >>  Tq ,  then  it  seems  reasonable  to  call  V(t) 
the  intensity  envelope,  or  merely  the  envelope.  If  Q(t)  is 
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II 

! 

I 


is  passed  through  a  lowpass  filter,  then  the  result  can  be  a 
good  estimate  of  the  sample  function  of  V(t).  As  stated  in 
the  introduction,  most  ambient  noise  experiments  concern  the 
squared  pressure  averaged  over  a  period  of  time  on  the  order 
of  a  minute.  The  resulting  time  series  can  be  sampled,  and 
the  samples  can  be  employed  to  obtain  estimates  of  the  statis¬ 
tical  parameters  of  the  envelope.  The  most  popular  statistic 
is  the  mean  value.  Greater  effort  is  required  to  estimate 
its  variance,  or  its  probability  distribution,  or  either  its 
covariance  or  power  spectrum. 

Estimates  of  the  amplitude  (not  intensity)  envelope  power 
spectrum  have  been  obtained  for  a  planar  array  with  720  direc¬ 
tional  elements  [10].  Two  one-octave  bands  were  processed: 

800  -  1600  Hz  and  1600  -  3200  Hz.  A  total  of  1500  beams  are 
DIKIJS  formed  to  cover  a  wide  sector  of  azimuth  and  elevation 
angles  between  plus  and  minus  90  degrees. 

The  amplitude  envelope  power  spectra  for  the  two  one-octave 
bands  are  shown  in  Figure  4.1  for  a  nearly  horizontal  beam. 

The  peak  above  30  mHz  is  attributed  to  an  oscillation  of  the 
array.  The  initial  negative  slope  is  about  8  dB  per  octave. 

It  is  considered  inadvisable  to  obtain  an  estimate  of  the 
autocovariance  function  by  Fourier  transforming  the  power 
spectrum.  However,  to  get  some  idea  of  the  temporal  implication 
of  this  data,  assume  a  first-order  Butterworth  spectrum 
(-6  dB  per  octave).  The  -3  dB  frequency  f  is  seen  to  be  about 
6  mHz,  and  the  associated  relaxation  time  is  d  *  (2  f0)_1-26  sec. 
The  relaxation  time  of  the  intensity  envelope  (amplitude  squared) 
is  half  of  that,  or  about  13  sec. 


LI 
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[GURE  4.1(a)  Low-band  Envelope  Spectrum-- 
Nearly  Horizontal  Elevation 
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FIGURE  4.1(b)  High-Band  Spectrum  of  a 

Nearly  Horizontal  Elevation 
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Figure  4.2  shows  the  amplitude  envelope  power  spectra 
in  the  two  one-octave  bands  for  a  beam  directed  at  the  sea 
surface.  Here  the  peak  in  the  vicinity  of  100  hZ  was  attri¬ 
buted  to  the  8-  to  10  second  period  wave  spectral  peak  of  the 
wind-generated  surface  waves.  Finally,  Figure  4.3  shows  the 
intensity  envelope  power  spectra  in  the  two  one-octave  bands 
for  the  output  of  a  single  (directional)  element.  Here  the 
energy  is  concentrated  at  lower  frequencies,  as  expected. 


For  the  multichannel  analog,  the  noise  input  may  be 
represented  by  a  vector  random  process 


Nj_(t)  =  /  V±(t )  Gi(t)  ,  i  =  1,2,3, - ,  n  (4.22) 


where  V^t)  is  a  non-negative  stationary  random  process 


with  average  value  . 


Gi(t)  is  a  zero-mean,  unit  variance  stationary  Gaussian 


process  statistically  independent  of  V^(t), 

j  =  1 .2 , 3, . ,  n. 


Now  define  ( t )  =  N^t);  then  the  joint  second  moment  function 


MQiJ (t)  u  E[Qi(t)  Qj (t  +  t  )] 


=  [ir^rrij  +  KVij(t)]  [1  +  ( 1 )]  (4.23) 


and  the  joint  covariance  is 


Kq i j ( t )  =  Ky^j(t)  +  2[m^mj  +  Ky^j(t)]  (4.^4) 
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[GURE  4.3(a)  Low-Band  Envelope  Spectrum- 
Single  Element 


4-12 


Report  No.  4262 


Bolt  Beranek  and  Newman  Inc. 


Assuming  that  TVi  >>  ,  all  i  and  all  j,  and  that 

Qi(t)  and  Qj(t)  are  subjected  to  lowpass  filtering,  then 

KVi j  (t  }  *  KQij(T)  (^-25) 


where  the  tilde  denotes  lowpass  filtering.  Nichols,  Young, 
and  Sayer  [11]  obtained  estimates  of  the  normalized  joint 
covariance  of  the  intensity  envelopes  of  pairs  of  beams;  i.e., 
of 


KVi,j(0) 


/  KVi(0)  KVj.  (0) 


(*.26) 


The  data  were  derived  from  a  bottomed  linear  array  operating  in 
a  six  Hz  band  dominated  by  shipping  noise.  Adjacent  beams 
overlap  at  -3  dB  response  levels.  In  the  azimuth  sector  with 
heavy  ship  traffic,  decreases  quite  rapidly  with  | i-j | , 

where  as  in  the  azimuth  sector  with  low  ship  traffic,  r^. 
decreases  rather  slowly  with  [ i- j | .  Because  of  their  security 
classification,  the  experimental  results  are  not  summarized 
herein. 


The  characteristics  of  received  signals  are  usually 
obtained  by  transmitting  signals  at  one  point  and  observing 
their  characteristics  at  another  point.  In  some  cases 
short  pulses  are  employed,  and  in  other  cases  CW  or  quasi  C\\' 
transmissions  are  employed.  The  latter  approach  Is  more 
suited  to  obtaining  the  envelope  of  a  narrowband  signal. 
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An  experiment  employing  both  types  of  signals  was 
conducted  in  water  with  an  average  depth  of  4  km. [12] 

Sources  and  receivers  were  suspended  at  a  depth  of  about 
1  km  below  the  tending  ships.  The  major  paths  were  refracted 
only;  the  upper  path  being  convex,  and  the  lower  one  being 
concave.  In  the  partially  saturated  regime  containing  the 
upper  path  the  amplitude  of  the  signal  received  from  a  quasi 
CW  transmission  was  expected  to  have  a  Rayleigh  distribution 
(or  an  exponential  distribution  of  intensity).  The  experimen¬ 
tal  data  were  consistent  with  that  expectation.  In  the 
unsaturated  regime  containing  the  lower  path,  the  amplitude 
was  expected  to  have  a  log-normal  probability  density  function 

f (a )  =  (2/  va)1/2  exp  [-2(lna  +  v/4)2]  (4.27) 

where  v  is  a  parameter  depending  on  the  frequency  of  the 
transmitted  signal.  In  the  frequency  region  from  1.76  to 
2.74  kHz,  the  parameter  is  given  by 

v  *  .35  f  (4.28) 

where  f  is  the  frequency  in  kHz. 

The  correlation  function  of  the  intensity  of  the  received 
signal  is  predicted  to  be  [13] 

<1 ( At )  I ( 0  )  >  =  1  +  |Q(At/to)|2  (4.29) 

where  the  average  intensity  is  <I>  =  1. 
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Q(  )  is  the  micropath  bandwidth  function.  Figure  4.4 
shows  the  predicted  correlation  function  (heavy  solid  line) 
and  measured  results  obtained  for  reciprocal  paths  (dashed 
and  solid  light  lines).  For  the  upper  path,  the  prediction 
and  the  measured  results  are  consistent,  and  the  decay  time 
is  about  5  minutes.  For  the  lower  path,  a  highly  stable 
reception  is  predicted  which  is  not  consistent  with  the 
measurements.  The  decay  time  of  the  latter  is  about  20  min. 

Measurements  spanning  a  period  of  about  a  half  year  were 
conducted  using  a  projector  operating  at  368  Hz  at  a  depth 
of  527  m  off  Eleuthera  and  a  pair  of  receivers  1318  km  away 
off  Bermuda  at  depths  of  1683  m  and  1723  m.[l4]  Average 
water  depth  between  the  locations  was  about  4500  m.  No  firm 
conclusion  was  reached  regarding  the  probability  density 
function  for  the  amplitude  envelope.  The  coefficient  of 
variation  (c^/m A)  for  the  amplitude  envelope  for  both  receivers 
averaged  0.51;  for  an  amplitude  envelope  with  a  Rayleigh 
distribution,  this  ratio  equals  0.523,  and  the  coefficient 
of  variation  for  the  intensity  (o^/m^)  equals  one.  Figure  4.5 
shows  the  normalized  autocovariance  function  for  amplitude 
fluctuations  for  one  of  the  receivers  derived  from  a  sixty-six 
hour  sample.  The  autocovariance  for  the  intensity  envelope 
is  the  square  of  that  for  the  amplitude  envelope.  For  the 
former,  the  1/e  decay  time  is  about  two  minutes. 

The  statistics  of  long-range  propagation  were  obtained 
by  having  sources  of  frequencies  of  10,  110,  and  262  Hz  at 
depths  of  110,  23,  and  25  m  respectively .[ 15 D  Signals  were 
recorded  from  both  bottom-mounted  and  suspended  receivers. 
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FIGURE  4.4  Intensity  correlation  for  the  upper  (top  figure) 
and  lower  (bottom  figure)  paths  at  2250  Hz. 


SOURCE  Ref  14 


FIGURE  4.5  Normalized  autocovariance  function  for  amplitude 
fluctuations. 
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Only  segments  with  a  mean  range  greater  than  600  km  and  a 
mean  range  rate  greater  than  10  km/hr  were  used.  Estimates 
of  Oj/m-j.  of  .89,  1.18,  and  1.14  were  obtained  for  the  fre¬ 
quencies  10,  110,  and  262  Hz  respectively.  These  ratios  are 
reasonably  consistent  with  a  Rayleigh  distribution  for  the 
amplitude  envelope  and  an  exponential  distribution  for  the 
intensity  envelope,  as  predicted  for  long  ranges  and  a 
range  spread  encompassing  at  least  one  convergence  zone. 
However,  neither  the  exponential  distribution  or  any  of 
several  other  candidate  distributions  produced  a  good  fit 
to  the  data  at  all  three  frequencies. 

Correlation  times  were  calculated  by  dividing  the 
duration  of  a  data  segment  by  the  number  of  independent 
samples  reported  in  [15].  Values  obtained  are  18.3,  4.5, 
and  1.9  min.  for  frequencies  of  10,  110,  and  262  Hz  respec¬ 
tively.  These  apply  to  range  rates  greater  than  10  km/hr. 

A  correlation  time  of  2  min.  was  reported  for  propagation 
at  368  Hz  from  a  fixed  projector  to  fixed  recei vers . [ 14 ] 

From  this  limited  base  of  data  it  can  be  concluded  that 
relaxation  times  of  both  signals  and  noise  can  range  from 
much  shorter  to  much  longer  than  a  typical  post -rect if icat ion 
averaging  time,  depending  on  the  circumstances. 
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5.0  FALSE  ALARM  PROBABILITY  -  COMMON  NOISE  POWER  ENVELOPE 

A  single  channel  (with  threshold  formation)  of  the  multi¬ 
channel  analog  is  shown  in  Figure  2.1,  and  an  equivalent  form 
is  shown  in  Figure  2.2.  If  the  noise  inputs  have  a  common 
power  envelope,  then  the  ith  sample  of  the  output  of  the 
summer  can  be  represented  by 

Yi  '  P1  tG0i  -  f  <GLn  +  Gn”  (5.1) 

-  PjUj  -  |Kj  )  (5.2) 

where  is  a  chi-squared  variate  with  one  degree  of 

freedom, 

is  a  chi-square  variate  with  two  degrees  of 
freedom,  independent  of  J^,  all  k. 

The  test  statistic  is  then 
n 


where  n  is  the  number  of  independent  samples. 

When  the  number  of  independent  samples  is  large,  then  the 
test  statistic  Z  is  Gaussian  given  {P^}  =  P,  and  the  approach 
discussed  in  the  Appendix  of  [1]  can  be  employed  to  obtain 
results  for  the  limiting  cases  of  fast  and  slow  envelope 
fluctuations.  This  approach  requires  the  mean  and  variance 
of  the  test  statistic  given  P.  These  can  be  readily  derived 
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using  the  recursion  formula  for  chi-square  variates: 

mn+1  =  (d  +  2n)mn>  (5.*0 

where  d  is  the  number  of  degrees  of  freedom. 

The  conditional  mean  of  (5-3)  is  then 

n 

E(Z | P)  =  (1-c)  £  P,  (5.5) 

1=1  1 

Squaring  (5*3)  gives 

Z2  =  (JiJj  -  cKiJi  +  )  ,  (5.6) 

and  taking  the  conditional  expected  value  gives 

E(Z2  |  P)  =  X)  P,2(3  -  2c  +  2c2) 
i  1 

+  ELPip/1  ‘  2c  +  c2)  (5.7) 

J 

-  X)  pi2(2  +  C2) 

+  EEPipi  (1  -  2c  +  c2  )  (5.8) 

i  J  J 

The  second  term  is  seen  to  be  the  square  of  the  conditional 

mean  as  given  by  (5.5);  thus  the  conditional  variance  is 

n 

Var  (ZjP)  =  (2  +  c2  )  £  P2  (5.9) 

i=l 
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The  conditional  mean  and  variance  can  be  expressed  more  compactly 
as 


E  (Z  | 

|P)  = 

(1  -  c  )  A 

(5.10) 

Var(Z 

|P)  = 

(2  +  c2  )B 

(5.11) 

where  A  =  ^  P^ 
i 

B  =  SPi 

i 

Then  (A-10)  of  [1]  can  be  expressed  as 


,oo  oo 


=  /da  / 
^0  •'0 


db  fAB(a,b)P 


Now  let  F  =  A/  /B  = 


{pi 


(l-c)a 
L  /(2+c2  )b  J 


(5.12) 


(5.13) 


then 


■/ 


df  fp(f)  P 


1-c 


/2  +  c' 


(5.1*0 


where  P(  )  is  the  normal  distribution  function. 


The  two  limiting  cases  can  be  readily  examined.  The 
first  is  the  limiting  case  of  very  slow  fluctuations  with 
*  p,  all  i.  Substitution  in  (5.13)  yields 

F  *  <=  /n  (5.15) 

/np’ 
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In  this  case,  the  density  function  for  F  is  a  delta  function 
at  F  =  /n;  evaluation  of  (5.14)  yields 


where  Q(  )  -  1  -  P()  =  P(-).  Since  the  power  envelope  is 
constant,  the  limiting  case  also  represents  one  with  Gaussian 
noise . 


The  second  is  the  limiting  case  of  very  fast  fluctuations. 
For  this  case  2Z  —  nm^. 


then 


n(mp + 


F 


/FT 


/n 


/l+(o  /m  ) 
P  P 


(5.17) 


where  op/mp  is  the  coefficient  of  variation  of  P.^ 

In  this  case  the  density  of  F  approaches  a  delta  function 
located  as  Indicated  by  (5-17).  Evaluating  (5.14)  gives 


=  Q 


( c- 1 )  /n 


•  /  ( 2  +  c  2  )(l  +  Op/rr^) 


(5.16) 
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If  the  power  envelope  is  a  chi-squared  process  then  (5.17)  can 
be  evaluated  by  (5.^): 


F  -+ 


d 

/(d+2  )d 


1 

/1+2/d 


and 


=  e 


( c  —  1 )/n 


/(2+c2 ) (1+2/d). 


(5.19) 


Table  5.1  shows  values  of  Pp  for  various  values  of  the 
threshold  parameter  c  and  the  number  of  degrees  of  freedom  d 
of  the  chi-squared  envelope.  The  number  of  averaged  indepen¬ 
dent  samples  is  10,000,  and  the  envelope  fluctuations  are  very 
fast.  At  low  threshold  settings,  the  variations  of  Pp  with 
the  degree  of  fluctuation  is  less  drastic  than  at  high  thres¬ 
hold  settings.  The  row  labelled  "  "  applies  to  both  Gaussian 
inputs,  and  very  slow  fluctuations  as  well.  As  expected,  the 
false  alarm  probability  increases  with  the  degree  of  envelope 
fluctuation,  which  varies  inversely  with  the  number  of  degrees 
of  freedom. 


TABLE  5.1 

False  Alarm  Probability  -  Fast  Fluctuations 


threshold  value  c 


d 

1.05 

1.07 

1.09 

2 

2.2  E-2 

2.6  E-3 

1.8  £-4 

4 

1.0  E-2 

6.3  E-4 

1.9  E-5 

8 

5.6  E-3 

2.1  E-4 

3.3  E-6 

m 

2.3  E-3 

3.9  E-5 

2.3  E-7 
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For  estimating  Pp  by  simulation,  the  analog  for  a 
channel  with  a  threshold  formed  from  two  neighboring  channels 
can  be  simplified  to  the  form  shown  in  Figure  5.1.  The  input 
labelled  represents  two  independent  squared  normal  inputs 
of  the  two  threshold  channels.  The  input  labelled  represents 
the  normal  input  to  the  center  channel. 

Table  5.2  shows  values  of  false  alarm  probabilities  for 
various  threshold  values  c,  and  various  noise  power  envelope 
relaxation  times.  The  number  of  terms  in  the  post-rectifica¬ 
tion  sum  is  60,  and  the  noise  power  envelope  is  a  chi-squared 
process  with  four  degrees  of  freedom.  For  the  limiting  case 
of  very  slow  fluctuations,  the  values  were  calculated  by  the 
method  developed  in  Section  3.0.  The  values  in  the  next  two 
rows  were  obtained  by  simulation  using  the  scheme  depicted 
in  Figure  5.1.  The  number  of  output  samples  was  10,000  in 
order  to  achieve  acceptable  accuracy.  Results  for  the  very 
fast  envelope  fluctuations  were  calculated  by  means  of  (5.19). 
In  all  columns  it  is  seen  that  the  false  alarm  probability 
varies  inversely  with  the  relaxation  time  of  the  noise  power 
envelope.  Also  note  that  the  variation  of  the  false  alarm 
values  with  relaxation  time  is  greater  for  larger  values  of 
the  threshold  value  c.  For  example,  the  ratio  of  Pp  for  fast 
fluctuations  to  that  for  slow  fluctuations  is  2.5  for  c  =  1.6, 
whereas  the  ratio  for  c  =  2.2  is  about  19-  It  is  the  latter 
case  that  is  more  relevant  to  sonar  operation. 
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TABLE  5.2 

False  Alarm  Probabilities  -  Variou s  Case s 


threshold  value  c 


envelope 

1.6 

1.8 

2.0 

2.2 

slow 

.0151 

.0033 

.0007 

.0001 

mN  *  60 

.0258 

.0072 

.0019 

.0003 

mN  '  20 

.032  6 

.0115 

.0041 

.0014 

fast 

.037  8 

.0135 

.0049 

.0019 

mN  is  the  number  of  samples  in  the  envelope 
N  formation  running  average. 


Section  4.2  discussed  beam  noise  measurements  obtained 
with  a  720-element  planar  array.  It  was  estimated  that  the 
relaxation  time  of  the  beam  noise  power  envelope  was  about 
13  seconds.  Relative  to  a  post-rectification  averaging  time 
of  two  minutes  or  120  seconds,  the  relaxation  time  is  indica¬ 
tive  of  a  fast  fluctuation.  Thus,  the  false  alarm  probability 
might  be  at  least  10  times  as  high  as  that  predicted  assuming 
either  very  slow  envelope  fluctuations  or  Gaussian  noise. 
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6.0  DETECTION  PROBABILITY  —  Common  Noise  Power  Envelope 

This  section  concerns  the  detection  of  signals  with  fluctuat¬ 
ing  power  envelopes  in  a  background  of  noise  inputs  with  a  common 
power  envelope.  The  schematic  of  the  simulation  for  this  case  is 
shown  in  Fig.  6.1.  The  schematic  for  the  generation  of  a  power 
envelope  is  shown  in  Fig.  2. A,  where,  for  this  case,  the  lowpass 
filters  LPF  are  running  averagers  RA. 

Results  give  the  probability  of  detection  as  a  function  of 
average  signal-to-noise  power  in  decibels.  Results  are  compared 
to  that  of  the  special  case  of  Gaussian  signals  and  noise. 

6.1  Gaussian  Signals  and  Noise 

To  provide  a  check  for  the  case  of  Gaussian  signals  and 
noise,  an  analytical  approximation  will  be  derived,  based  on 
the  assumption  that  the  test  statistic  is  nearly  Gaussian.  In 
that  case,  (2.13)  of  Ref.  1  becomes 


P 


D 


dx  Z(x) 


(6.1) 


where  Czs  is  the  coefficient  of  variation  of  the  test 
statistic  when  the  signal  is  present 


Z (  )  is  the  Gaussian  (0,  1)  density  function 

The  inverse  of  C^s  is  given  by  (6.5)  of  Ref.  1;  for  the  sampled 
data  case,  replace  the  bandwidth-time  product  WT  by  n/2 ,  where  n 
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is  the  number  of  independent  samples  added  by  the  post¬ 
rectification  running  summer.  And  for  a  two-channel  threshold, 
let  ci  =  c/2;  then  (6.5)  of  Ref.  1  becomes 


r  +  1  -  c 


r  +  l)2 


+  c2/2 


(6.2) 


where  r  =  p  /p  is  the  signal-to-noise  power  ratio.  Substituting 
(6.2)  in  (6.1)  gives 


pd  -  Q 


N  n/2  (c  -  1  -  r) 


^  2  (r  +  1)  +  c2/2 


(6.3) 


Fig.  6.2  shows  the  transition  curve  for  n  *  60  calculated 
from  (6-3).  The  detection  differential  (10  log]_o  r  for  Pd  =  0,5) 
is  .8  dB,  and  the  spread  (dB  between  PD  =  0.2  and  PD  =  0.6)  is 
3.05  dB.  For  a  single-channel  analog  with  a  "linear"  detector, 
a  deterministic  threshold,  and  a  sine  wave  signal  in  stationary 
Gaussian  noise  [16],  the  detection  differential  is  -2.7  dB,  ana 
the  spread  is  2.4  dB  for  n  =  64.  For  the  same  type  of  detector 
and  inputs,  the  detection  differential  is  -13-7  dB  and  the  spread 
is  2.05  dB  for  n  =  8192. 


It  is  clear  from  the  horizontal  axis  of  Fig.  6.2  that  the 
small-signal  regime  does  not  obtain  for  n  =  60,  with  the  consequence 
that  the  transition  curves  that  follow  cannot  be  extrapolated  to 
larger  values  of  n.  The  simulation  analog  was  run  with  Gaussian 
signals  and  noise  for  n  =  60,  and  the  number  of  output  samples 
was  1000.  The  results  Ps  of  the  simulation  are  compared  with 
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the  results  PA  obtained  using  the  approximation  (6.3)*  The 
differences  are  seen  to  be  small,  and  the  consistent  sign  of 
the  differences  indicates  a  bias  in  the  approximations  calculated 
from  (6.3). 

Table  6.1 


Comparison  of  Simulation  and  Calculation 


10  1 oq  r ,  dB 

Is  ' 

Ia 

PS  '  PA 

-3 

.039 

.038 

+  .001 

-2 

.089 

.083 

+  .006 

-1 

.188 

.175 

+  .013 

0 

.346 

.333 

+  .013 

1 

.558 

.  547 

+  .011 

2 

.776 

.756 

+  .018 

3 

.916 

.902 

+  .014 

4 

.980 

.969 

+  .011 

6.2  Limiting  Case  —  Very  Fast  Fluctuations 

This  section  derives  an  analytical  result  for  PD  for  the 
limiting  case  of  signals  and  noise  with  fast-fluctuating  power 
envelopes.  The  procedure  followed  is  the  same  as  that  employed 
in  Section  5-0  to  obtain  P^.  The  procedure  requires  the  mean 
and  variance  conditioned  on  random  variates  derived  from  the 
power  envelopes.  The  conditional  means  and  variances  will  be 
derived  for  the  continuous  input  case  and  then  converted  to  the 
sampled  data  case. 

Equation  ( A— 3 )  of  Ref.  1  gives  the  conditional  mean  for  a 
general  case  with  continuous  inputs.  For  the  case  of  a  common 
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noise  power  envelope  and  a  two-channel  threshold  with  weighting 
factor  c^  =  c/2,  the  conditional  mean  is 

E[Z | P(t )  ]  -  As  +  (1  -  c)An  (6.4) 

where  _P(t)  is  the  set  of  power  envelopes 

r 

T_1  I  du  Pg(u) 

0 

post-rectification  averaging  time 

AN  =  T_1/  du  PN(u) 

0 

Equation  (A-9)  of  Ref.  1  gives  the  conditional  variance  for  a 
general  case  with  continuous  inputs.  For  the  specific  case 
described  above,  the  conditional  variance  is 


AS 


T  is 


Var  [2 | P(t ) ]  =  (WT)'1  [Bg  +  2BC  +  (1  +  c2/2)BN]  (6.b) 


where  W  is  the  input  bandwidth 


Bo  = 


B„  = 


T-i  j 1  du  P|(u) 
0 


•/ 


T'1  /  du  Ps(u)PN(u) 


B. 


*  T”  lj" 


du  P2(u) 
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For  the  case  of  fast  fluctuations,  examination  shows  that 
-*■  ps,  Bg  p|  +  o|,  etc.;  and  if  the  WT  product  is  sufficient 
the  distribution  approaches  normal.  For  the  continuous  input 
case,  the  detection  probability  is  approximately 


>JW[(c  -  1)  pN  -  Pg] 

PS  +  °S  +  2pSPN  +  (pN  +  (1  +  cV2) 


(6.6) 


For  the  sampled-data  case,  the  number  of  independent  samples 
that  are  averaged  is  n  =  2WT;  thus 


VQ 


A 


(c  -  1  -  r) 


^  r2  (1  +  C2)  +  2r  +  (1  +  Cj)  (1  +  c2/2) 


(6.7) 


where  r  =  p0/p 
s  N 


cs  =  °s/ps 


CN  *  VPU 


Finally,  if  the  envelopes  are  chi-squared  processes  then 
C2  =  2/dg,  C2  =  2/d.,,  and 


f  (c  -  1  -  r) 

1 

_ 

IJr2  (1  +  2/d,,)  +  2r  +  (1  +  2/d,' 

(1  +  c2/2 ) 
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Comparison  of  (6.7)  with  (6.3)  for  the  Gaussian  input  case 
shows  that  the  denominator  of  the  argument  is  larger;,  hence,  the 
transition  curve  for  the  fast  fluctuation  case  will  have  a  larger 
spread  in  signal-to-noise  power  ratio. 


6.3  Limiting  Case  -  Very  Slow  Fluctuations 


This  section  derives  an  analytical  result  for  for  the 
limiting  case  of  signals  and  noise  with  signal  and  noise  power 
envelopes  that  are  very  slowly  fluctuating.  In  the  limit,  the 
power  envelopes  are  random  variates  vice  random  processes.  For 
this  case,  the  general  result  is* 


-f 

0 


dr  fR(r)  g(r ) 


where  fR(  )  is  the  density  of  R 


(6.9) 


R  =  Ps/  Rv; 

g(  )  is  the  transition  curve  for  Gaussian  inputs. 


If  the  envelope  variates  are  gamma-distributed,  then  the 
probability  density  of  R  is  given  by  (6.27)  of  Ref.  1.  For 
the  particular  case  of  chi-squared  variates,  then  the  prob¬ 
ability  density  is 


where  r  =  ps/pN 
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Given  an  adequate  time-bandwidth  product,  then  the  transition 
curve  for  Gaussian  inputs  is  given  by  (6.3).  Thus,  for  the 
limiting  case  of  very  slowly  fluctuating  chi-squared  power 
envelopes,  the  probability  of  detection  is  found  by  substituting 
(6.3)  and  (6.10)  in  (6.9). 

6.4  Transition  Curves 

Transition  curves  obtained  by  both  calculation  and  simula¬ 
tion  are  shown  in  Fig.  6.3*  The  transition  curve  for  Gaussian 
inputs,  discussed  in  Section  6.1,  shows  the  least  spread  in  sigr.ai- 
to-noise  power  ratio.  All  the  other  curves  pertain  to  chi- 
squared  envelopes  with  two  degrees  of  freedom  for  the  signal 
(appropriate  for  propagation  over  a  large  range),  and  four 
degrees  of  freedom  for  the  noise.  These  curves  clearly  show 
that  the  faster  the  fluctuation,  the  smaller  the  spread  in 
signal-to-noise  power  ratio.  The  curves  diverge  more  at  the 
positive  signal-to-noise  levels  than  at  the  negative.  At  the 
positive  levels,  the  signal  power  dominates,  and  in  this  case, 
it  has  a  greater  variability  than  the  noise. 

The  solid  and  dashed  lines  in  Fig.  6.4  repeat  two  of  the 

transition  curves  shown  in  Fig.  6.3:  n^  =  20,  n„  =  20  and  ng  = 

60,  n^  =  60.  The  triangles  mark  points  obtained  by  simulation 

for  the  case  n^  =  20,  n^  =  60.  For  R  ^  1  dB,  these  points  fail 

very  near  to  the  n^  =  20,  n^  =  20  curve;  this  is  reasonable 

since  both  cases  have  the  same  signal  power  envelope  relaxation 

time.  For  R  <  -  1  dB,  these  points  fall  nearer  the  ng  =  60, 

n,,  *  60  curve;  this  is  reasonable  since  these  cases  have  the 
N 

same  noise  power  envelope  relaxation  time.  The  circles  mark 
points  obtained  by  simulation  for  the  case  n^  =  60,  n^  =  20. 

For  R  >  2  dB,  these  points  fall  very  close  to  the  curve  for 
the  ns  =  60,  nN  *  60  case.  For  H  <_  -2  dB,  these  points  fall 
very  near  the  curve  for  the  ns  ■  20,  n^T  =  20  case. 


6-9 


Newman  Inc 


in 


i 


cc 

=> 

CD 


I 


i 


Report  No.  4262 


Bolt  Beranek  and  Newman  Inc. 


r 


Each  transition  curve  can  be  partially  characterized  by 
a  pair  of  numbers:  R^,,  the  value  of  R  for  which  =  0.5; 
and  the  spread  in  dB  between  R  values  corresponding  to  P^  = 

0.2  and  P^  =  0.8.  Threshold  and  spread  values  are  shown  in 
Table  6.2,  which  shows  that  the  range  of  R^  is  only  0.7  dB, 
whereas  the  range  of  spreads  is  almost  8  dB. 

Historically,  predictions  of  detection  performance  usually 
fall  into  one  of  two  categories.  One  of  them  assumes  that  the 
inputs  are  Gaussian,  and  the  other  assumes  that  the  relaxation 
time  of  the  inputs  is  very  large  compared  to  the  post-reetif ieacic 
averaging  time.  These  assumptions  lead  to  the  extreme  transition 
curves  labeled  G  and  S  on  Fig.  6.3  respectively.  The  latter 
assumption  could  lead  to  excessive  optimism  regarding  the  detec¬ 
tion  of  sources  at  long  ranges  from  multiple  independent  obser¬ 
vations.  The  former  assumption  could  lead  to  excessive  optimism 
regarding  detection  performance  at  short  ranges. 


Table  6.2 


Threshold 

and 

Spread 

Values  - 

Common 

No  i ; 

se  Power 

En  ve  1  ope 

Si  qnal 

Noise 

*T 

,  dB 

Spreac 

i,  dE 

Gaussian 

Gaussian 

0 

.8 

3. 

1 

*  2,  fast 

ds 

*  4  ,  fast 

0 

.6 

3. 

.  7 

'  •  "S  * 

20 

"  n  « 

.  nN 

20 

0 

.  6 

6. 

,  3 

"  •  ns  ' 

20 

"  n  * 

.  nN 

60 

0 

.6 

5. 

,C 

.  ns  * 

60 

*  nN  * 

20 

1 

.  3 

6. 

,  3 

‘  •  "S  * 

60 

•  nS  * 

60 

1 

.  1 

7, 

,2 

"  .slow 

"  ,  slow 

1 

.5 

11. 

,0 

!| 
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7.0  PARTIALLY-CORRELATED  NOISE  POWER  ENVELOPES 
7.1  Introduction 

The  previous  section  dealt  with  a  case  in  which  a  noise 
power  envelope  was  common  to  the  inputs  of  all  of  the  channels. 
This  is  equivalent  to  saying  that  the  noise  power  envelopes 
are  fully  correlated.  This  section  deals  with  a  case  in 
which  the  noise  power  envelopes  are  partially  correlated. 

Ref.  11  presents  data  showing  the  correlation  of  averaged 
power  levels  of  beam  outputs  derived  from  a  line  array  operat¬ 
ing  in  a  frequency  band  in  which  the  dominant  noise  is  gener¬ 
ated  by  ships.  The  data  are  not  repeated  here  because  of 
their  security  classification. 

Figure  7-1  illustrates  a  means  for  generating  three 
partially-correlated  chi-squared  envelopes  with  four  degrees 
of  freedom.  The  correlation  between  the  center  channel 
envelope  and  that  of  a  threshold  channel  is  0.75,  and  the 
correlation  between  the  threshold  channels  is  0.5.  For 
greater  generality,  the  running  averagers  (RA)  could  be 
replaced  by  linear  lowpass  filters. 

Figure  7.2  illustrates  the  flow  diagram  for  the  multi¬ 
channel  analog  and  its  input  components  appropriate  for 
obtaining  transition  curve  data.  False  alarm  probabilities 
can  also  be  estimated  by  setting  the  parameter  k  equal  to 


zero. 


r 
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FIGURE  7.1  Synthesis  of  Partially-Correlated  Noise 
Power  Envelopes. 
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7.2  Limiting  Case-Very  Slow  Fluctuations 

The  limiting  case  of  very  slow  fluctuations  will  be 
examined  first  with  the  objective  of  predicting  false  alarm 
probability.  In  this  case,  the  power  envelopes  are  random 
variates  rather  than  random  processes.  A  single  channel 
(with  threshold  formation)  of  the  multichannel  analog  is 
shown  in  Figure  2.1,  and  an  equivalent  form  is  shown  in 
Figure  2.2.  With  power  envelopes  described  above,  the  test 
statistic  can  be  expressed  as 


n 


z  ■ 

2  (PlGli  + 

p-iG-n):i 

(7.1) 

i=l 

n 

n 

n 

-  po£°oi  - 

2  (P1  ^  Gli 

+  p.a  E  Gijj) 

1  i=l  21 

(7.2) 

where  the  subscript  0  denotes  the  center  channel.  The  sums 
shown  in  (7.2)  represent  chi-square  variates  of  n  degrees 
of  freedom.  As  n  increases,  the  distribution  of  the  chi- 
squared  variate  approaches  normal.  If  the  sums  are  assumed 
normal,  then 


Z  =  PnGn  -  £  ( P,  G ,  +  P  ,  G  ,  ) 


0  0 


ri 


-i-i 


(7.3) 


where  G^  is  a  Gaussian  random  variate  with  mean  n 
and  variance  2n,  i  =  -1,  0,  1 


If  the  envelopes  are  syntehsized  in  a  manner  analogous 
to  Figure  7.1  then 


5  ^  6 

Z  a  GnEKi  '  f  (Gi  S  +  G  -,  £  )  (7.4) 

0  i=2  1  6  1  i=l  i=3 
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where  Ki  = 

GEi  =  Gauss  (0*1) 

Finally,  Pp  =  P(Z  >  0)  (7.5) 

Equation  (7.4)  and  (7-5)  form  the  basis  for  estimating  Pp 
by  means  of  a  Monte  Carlo  simulation. 

Figure  7.3  shows  the  false  alarm  probability  for  a  single 
observation  as  a  function  of  the  threshold  parameter  c  for 
both  partially-  and  fully-correlated  slowly-fluctuating  noise 
power  envelopes.  The  difference  of  the  false  alarm  probabilities 
is  seen  to  increase  with  threshold  value.  At  c  =  2.2,  the 
partially  correlated  noise  power  envelope  yields  a  false 
probability  that  is  almost  100  times  that  produced  by  the 
fully  correlated  envelope. 

7.3  Limiting  Case  -  Very  Fast  Fluctuations 

A  means  will  be  derived  for  calculating  the  false  alarm 
probability  for  the  case  of  partially-correlated  noise  power 
envelopes  with  very  fast  fluctuations.  The  approach  discussed 
in  the  Appendix  of  [1]  will  be  employed.  For  the  work  at  hand, 
(A-10)  of  [1]  can  be  written  as 

Pp  =  dr  fR(r)  P (r )  (7.6) 

where  P()  is  the  distribution  function  for  a  normal  (0,1) 
variate , 
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fR()  is  the  probability  density  function  for  R, 

R  =  E (Z } P) 

Vv(  Z  IP) 

V()  denotes  the  variance  of  the  argument 

P  =  { }  denotes  a  set  of  sample  values  of  the 

noise  power  envelopes. 


For  this  case,  the  test  statistic  can  be  expressed  as 
n 

z  -  I>0iGoi  -  I  <pnGh  +  ”-e> 


The  expected  value  of  the  test  statistic  Z  given  the  values 
denoted  by  P  is 

n 

E(z  v  -  E  cp0i  -  I  <pn  +  p-n>:>  (7-9) 

If  the  envelopes  are  synthesized  according  to  the  scheme 
depicted  by  Figure  7.1,  and  if  only  one  term  is  employed  in 
the  running  averagers  (this  corresponds  to  the  most  rapid 
fluctuation),  then  the  jth  sample  of  the  power  envelope  for 
channel  i  can  be  expressed  as 


(7.10) 


where  is  the  jth  sample  of  input  k,  a  normal  random 

variate  (0,1),  statistically  independent  of  H, 


unless  k  =  p  and  j  =  q. 


pq 
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Substituting  (7.10)  in  (7*9)  and  collecting  terms  gives 


n 

E(ZjP)  =  £“|  Hli  +  (1"l)H2i  +  (1~c)H3i 


i=l 


+  (l-c)Hji  +  (l-f)H^  -  §-H|a] 


(7.11) 


If  n  is  large,  then 
n 

-  n,  j  =  1,2, 3,  ^>5, 6 

i=l  J 

and  E ( Z | P )  *  4n(l-c) 


(7.12) 


The  conditional  variance  of  the  test  statistic  is  found 
by  obtaini:  -  the  conditional  mean  square  and  subtracting  the 
square  of  the  conditional  mean.  The  result  obtained  is 

n  2 

v(z]p)  =  2  23cP20i  +  $  (P'n  +  Pin)]  (7.13) 

Squaring  and  summing  over  n  samples  gives 


i+5  i+5  n 


-2 

k=i+2  £=i+2  3=1 


2  H2 


(7.1*0 


If  n 


is  large. 


then 


i  +  5  i  +  5 

i  +  5 

=  n 

k= 

£  5U 

k=i+2 

k^S. 

i  +  5  i+5 

i+5 

=  n 
k 

EEi 

=i+2  fi.=i+2 

k=i+2 

=  24n 


(7.15) 


(7.16) 
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(7.16) 


Substituting  (7.12)  and  (7.18)  in  (7.7)  and  simplifying  give: 


R  _  /n(l-c ) 

/(2+ci ) 3/2 


(7.19) 


The  density  function  for  R  is  approximately  a  delta  function 
at  the  location  indicated  by  (7.19).  Evaluation  of  (7-6)  yields 


«  Q 


/n(l-c) 
'(2+c*  )  3/2 


(1.20) 


The  result  is  identical  to  (5-19)  with  dJv,  *  4.  That  is,  for 
fast  fluctuations,  the  false  alarm  probability  is  independent 
of  the  degree  of  envelope  correlation. 

7.4  False  Alarm  Probability 

Table  7.1  shows  values  of  flase  alarm  probabilities  for 
various  threshold  values  c,  and  various  noise  power  envelope 
relaxation  times.  The  number  of  terms  in  the  post-rectification 
sum  is  60,  and  each  noise  power  envelope  is  a  chi-squared 
process  with  four  degrees  of  freedom.  For  the  limiting  case 
of  very  slow  fluctuations,  the  values  were  generated  by  the 
Monte  Carlo  method  described  in  Section  (7.2).  The  values  in 
the  next  two  rows  are  obtained  by  simulation  as  described  in 
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TABLE  7.1 

False  Alarm  Probability  -  Partially  Correlated  Envelope 


value  of  c.  threshold  parameter 


envelope 

1.6 

1.8 

2.0 

2.2 

slow 

.0733 

.03  65 

.0212 

.0123 

nN  *  60 

.0674 

.0305 

.0141 

.0063 

mN  *  20 

.0510 

.0203 

.0078 

.0046 

fast 

.0378 

.0135 

.0049 

.0019 

Section  7.1.  The  number  of  output  samples  was  10,000  in  order 
to  achieve  acceptable  accuracy.  Results  for  the  very  fast 
envelope  fluctuations  were  calculated  by  means  of  (7.20).  In 
all  columns  it  is  seen  that  the  false  alarm  probability  varies 
directly  with  the  relaxation  time  of  the  power  envelope.  This 
is  just  the  reverse  of  the  case  of  the  fully-correlated  noise 
power  envelopes  (i.e.,  the  common  noise  power  envelope). 
Comparison  of  the  last  rows  of  Tables  7.1  and  5.2  show  common 
values  of  false  alarm  probabilities  for  the  limiting  case  of 
fast  fluctuations  for  the  two  classes  of  envelopes.  Again, 
it  is  noted  that  the  variation  of  false  alarm  probability  with 
relaxation  time  is  greater  for  larger  values  of  the  threshold 
value  c.  For  example,  the  ratio  of  Pp  for  slow  fluctuations 
to  that  for  fast  fluctuations  is  1.9**  for  c  =  1.6,  whereas  the 
ratio  for  c  =  2.2  is  about  6.5* 
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The  greatest  disparity  between  the  values  in  Tables  7.1 
and  5*2  occurs  for  the  limiting  case  of  very  slow  fluctuations. 
Recall  also  that  the  first  line  of  Table  5-2  also  corresponds 
to  a  Gaussian  noise  input.  For  a  threshold  value  of  2.2,  the 
false  alarm  probability  for  the  partially  correlated  noise 
power  envelope  is  about  120  times  that  for  a  fully-correlated 
noise  power  envelope  or  for  Gaussian  noise  inputs. 

7.5  Transition  Curves 

Four  transition  curves  obtained  by  simulation  are 
represented  in  Figure  7*4.  All  curves  pertain  to  chi-squared 
envelopes  with  two  degrees  of  freedom  for  the  signal  (appropriate 
for  propagation  over  a  long  range),  and  four  degrees  of  freedom 
for  the  noise.  The  dashed  and  solid  curves  pertain  to  ng=20, 
n^  =  20,  and  ng  =  60,  n^  =  60  respectively.  Of  these,  the 
former  shows  the  least  spread.  The  triangles  mark  points 
obtained  for  the  case  ng  =  20,  =  60.  For  R  >_  0  db,  these 

points  fall  near  the  ng  =  20,  nf.  =  20  curve;  this  is  reasonable 
since  both  cases  have  the  same  signal  power  envelope  relaxation 
time.  For  R  >_  -2  dB,  the  points  in  the  triangles  approach  the 
ng  -  60,  n^  =  20  curve;  this  is  expected  since  these  cases 
have  the  same  noise  power  envelope  relaxation  time.  Similar 
remarks  can  be  made  about  the  cases  represented  by  the  solid 
curve  and  the  circled  points.  In  these  respects  the  transi¬ 
tion  curves  for  the  cases  of  partially-correlated  noise  power 
envelopes  are  similar  to  those  for  the  common  noise  power 
envelope . 


7-n 
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Table  7-2  gives  values  for  R^,  the  value  of  R  for  which 
PD  =  0.5,  and  the  spread  in  dB  between  R  values  corresponding 
to  =  0.2  and  P^  =  0.6. 


TABLE  7.2 

Threshold  and  Spread  Values-Partially-Correlated 
Noise  Power  Envelopes 

Signal,  ds«2  Noise,  d$«4  dB  spread,  dB 


20 

20 

0.0 

5.7 

20 

60 

-0.2 

6.2 

60 

20 

0.7 

6.8 

60 

60 

0.4 

7.2 

The  range  of  threshold  values  (maximum  minus  minimum  value) 
is  only  0.6  dB,  whereas  the  range  of  spread  values  is  1.5  dB. 
Comparing  the  entries  in  Table  7.2  to  the  corresponding  entrie 
in  Table  6.2  will  not  be  highly  informative  since  the  values 
of  the  threshold  parameters  were  different  (2.0  vice  2.2, 
respectively).  Qualitatively,  the  transition  curves  for  the 
two  cases  (partially  and  fully  correlated  noise  power  envelope 
are  similar. 
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APPENDIX  A 

Moments  of  the  Test  Statistic  and  Pp  Predictions 


An  analytical  method  was  devised  to  calculate  the  detection 
performance  of  the  multichannel  analog  operating  with  fluctuating 
inputs.  A  means  was  found  to  approximate  the  density  function  of 
the  test  statistic  by  an  expansion  of  the  form 

n  i 

f7  (x)  =  g(x)  I  a  .  x  (A-l) 

L  i  =  l 

where  g(x)  is  a  basis  density  function.  The  coefficients  ai  of 
the  expansion  are  determined  by  requiring  that  its  first  n  moments 
match  those  of  the  test  statistic.  The  resulting  set  of  n  +  1 
equation  is 

if0  aink+i  =  ak’  k  =  ^  1»  n  (A'2) 

where  nu  is  the  i£il  moment  of  the  basis  density  function 
a,  is  the  k—  moment  of  the  test  statistic 

■K 

Moments  to  the  third  order  were  previously  derived  [1] 
and  applied  to  predict  detection  performance.  More  recently, 
the  fourth-order  moment  was  derived,  and  the  moments  were  applied 
to  predict  false  alarm  probability.  Results  of  those  calculations 
are  presented  in  the  sequel,  along  with  the  fourth-order  moment 
of  the  test  statistic,  shown  below. 


A- 1 
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y4z  = 


■where 


12(pk6)1'  ♦  Us-1  Iv)  u 

+  S-1  IDS  [ 2n-1  ID,,  +  (l+2n‘1)y”2]  u2 


+  n  (idi;  + 


1  . -2-1 


+  ^2idn  +  l6n_1  hn  +  ^  JDN  <1+2W  +  l6rr"(2Iih  +  V] 

+  ^Y^  [l+fe^d+PI^,)  +  32n'2  IDN  +  4n"2  d+21^,)  +  l6n'3  (21^.41^0} 


-1 


-2/ 


(.  A  —  o  ) 


PN  is  average  noise  power 


y  is  the  false  alarm  parameter.  For  a  normal  input,  the 
probability  of  false  alarm  is  P(-y)  *  Q(y). 

P(  )  is  the  cumulative  distribution  for  a  zero-mean, 
unit  variance  normal  random  variable, 

Q(  )  *  1  -  P(  ) 

m  is  the  number  of  channels  employed  for  threshold 
formation, 

W  is  the  input  bandwidth  of  the  receiver, 

T  is  post-rectification  averaging  time, 

s  is  the  number  of  degrees  of  freedom  of  the  signal 
power  envelope. 

the  subscripts  S  and  N  denote  signal  and  noise 

respectively, 

1^  «  T-2y* du J dv  r2(u-v)  (A-4) 

I,p  ■  T-3  f du  /*dv f  dw  r(u-v)r(v-w)r(w-u)  (A-3) 


A-2 
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I"1‘/du/dv/d“/  dx 


r(u-v)r (v-w)r (w-x )r (x-u  ) 


(A -6) 


the  limits  of  integration  are  from  0  to  T. 

r()  is  the  autocovariance  function  for  the  (Gaussian) 
components  of  the  chi-squared  envelope  processes, 
in  is  the  autocovariance  function  for  the  (Gaussian) 
components  of  the  chi-squared  envelope  process. 

The  integrals  defined  in  the  previous  paragraph  were 
derived  for  the  case  in  which 


r (t )  =  exp  t-|t|  (2D)"*1] 

where  D  is  the  relaxation  time  of  the  power  envelope.  The 
results  are 


2t  t1  + 

?  (e'T/D  -  1)] 

U-8) 

-2§  +  (1  +  2^)e~T/D] 

(A-9) 

12(jf)2  \ 

[2|  -  6(^)2  +  [1  +  4^  4  6(|)2]e_T/D} 

( A-1Q ) 

The  results  were  initially  applied  to  predict  the  false 
alarm  probability.  For  the  case  in  which  the  noise  power 
envelope  had  a  low  coefficient  of  variation,  the  normal 
density  was  employed  for  the  basis  function,  and  the  false 
alarm  probability  is  then  given  by 
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where 


a4 


is  the  skewness  coefficient  of  the  analog  output. 


is  the  kurtosis  of  the  output,  i.e.,  t 


Table  A-l  shows  results  for  a  case  in  which  the  noise  power 
envelope  is  a  chi-square  process  with  100  degrees  of  freedom. 
The  table  also  shows  P*.  for  a  Gaussian  input  in  the  second 

r 

column . 

Table  A-l.  False  Alarm  Probability 


False  alarm  probability 
Gauss.  Chi-squared  envelope,  d.f.  *  100 


Y 

— 

D/T  ■*« 

r— 

1! 

H 

Q 

D/T  =.3 

D/T  =  .1 

D/T  =.03 

D/T  =0 

0.5 

.3085 

.3085 

.3089 

.3095 

.3100 

.3102 

.3103 

1.0 

.1587 

.1585 

.1592 

.1600 

.1606 

.1609 

.1611 

1.5 

.06681 

.06681 

.06730 

.0679 

.0684 

.0636 

.0687 

2.0 

.02275 

.02291 

.02312 

.0234 

.0236 

.0238 

.0235 

2.5 

.006210 

.00633 

.00640 

.00649 

.00658 

.00663 

.00665 

3.0 

.001350 

.001355 

.00139 

.00143 

.00146 

.00148 

.00149 

For  the  case  D/T-*-®,  the  power  envelope  is  a  random  variable 
rather  than  a  random  process,  and  it  can  be  shown  that  P?  is 
Independent  of  the  envelope  distribution,  and  that  it  is 
equal  to  that  for  a  Gaussian  input.  The  Table  shows  good 
agreement  between  the  values  in  the  Gauss  and  D/T+®  columns. 
Since  the  difference  in  output  moments  of  second  or  higher 
order  are  the  most  pronounced  for  these  two  cases,  it  is  con¬ 
cluded  that  the  figures  in  the  other  columns  are  also  accurate. 


A-4 
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For  values  of  y  greater  than  3-5,  however,  the  predicted  value 
of  Pp  for  D/T-*-®  diverge  appreciably  from  the  Gauss  input  values. 
This  can  be  anticipated  from  the  small  number  of  terms  in  the 
expansion  and  the  distance  from  the  center  of  the  distribution. 
The  remaining  columns  show  that  Pp  increases  inversely  with  the 
D/T  ratio. 

The  right-most  column  of  Table  A-l  shows  Pp  for  the  case  of 
very  fast  power  envelope  fluctuations.  These  results  were  cal¬ 
culated  from  (5.19). 

The  case  of  the  low-fluctuation  envelope  is  important  from 
the  theoretical  point  of  view,  but  of  limited  interest  from 
the  operational  point  of  view  because  of  the  very  small  changes 
of  probability  values. 

For  more  general  application,  the  basis  density  employed 
is  that  for  the  output  (denoted  by  PW)  of  the  analog  when  the 
power  envelope  is  a  random  variable.  For  this  case,  the 
expression  for  the  false  alarm  probability  is 


where 


£  a± 

i-1 


miP  Ii 


(A-12) 


a.j.  is  the  i  weighting  coefficient  found  by  solving 
(3.6)  of  Ref.  1. 

mip  is  the  ith  moment  of  the  envelope  variate 

f  ()  is  the  normal  density  function  with  ,-ean 

w  _  _2 

m  *  -B  Vl+n-1,  and  variance  - 


1  w  J-0  J  J 
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A  recursion  formula  was  derived  for  Lj  : 

Lj  =  yj_1  Z(y)  +  (J-1)Lj_2  (A-13) 

where  Z(y)  is  the  normal  density  function  for  zero  mean,  unit 
variance 


Lo  =  Q(y)  (A-U) 

The  results  were  applied  to  a  case  in  which  the  power  envelope 
process  is  chi-squared  with  100  degrees  of  freedom,  and  y  =  1- 
The  results  are  shown  in  Table  A-2,  which  also  shows  results 
using  the  normal  basis  for  the  same  case. 


Basl  s 
PW 

Gauss 


Table  A-2 

PF  calculations  for  n*100,  y*1 
False  alarm  probability 

Gauss  chi-squared  envelope,  d.f.  ■  100 _ 

D/T  D/T  =  l  D / T=1  / 3  D/T*.3  0/T=.l  D/T*  0 

.1563  .15865  .1587  .1591  -  .1588  (.1611) 

.15865  .1585  .  1592  . 1600  .1606  (.161  1  ) 


Examination  of  Table  A-2  shows  that  the  gaussian  basis  yields  more 
accurate  results  than  the  PW  basis. 


The  PW  basis  was  also  applied  to  cases  of  high  variability 
and  was  found  to  be  unsuitable.  Therefore,  consideration  was 
given  to  the  use  of  a  simulation  or  Monte  Carlo  approach 
based  on  the  multichannel  analog  and  the  compound  random, 
process.  This  approach  is  covered  in  the  text  of  this  report. 
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APPENDIX  B 


Multi-Channel  Simulation  Software 


The  software  for  the  multi-channel  simulator  is  composed  of 
seven  Fortran  programs.  These  programs  can  be  separated  into 
two  groups:  random  number  generators  and  processor  models. 

The  random  number  generators  generate  the  random  processes  used 
as  input  for  the  processors.  The  processes  series  are  written 
to  disk  as  an  intermediate  step  so  that  they  may  be  used  mere 
than  once  by  the  processor  programs. 

B-l.  Random  Number  Generator  Programs 

There  are  three  random  number  generator  programs,  corres¬ 
ponding  to  Gaussian,  exponential  and  Chi-squared  random  variables. 
Each  of  the  programs  receives  generation  instructions  from  input 
at  logical  I/O  unit  5.  Each  program  writes  out  one  file  of  random 
variables  for  every  two  lines  of  input  from  unit  5.  The  first 
line  assigns  values  to  variables  in  the  Fortran  namelist  called 
FACT  (see  a  Fortran  manual  for  an  explanation  of  the  namelist 
I/O  feature).  This  line  must  begin  with  a  blank  immediately 
followed  by  "$fact."  The  line  must  end  with  a  with  the  input 
variable  assignments  falling  in  between.  (Actually,  the  name¬ 
list  input  may  extend  over  an  arbitrary  number  of  lines,  but  must 
begin  and  end  as  described.)  After  namelist  input,  the  program 
reads  the  name  it  will  assign  to  the  file  it  is  about  to  generate. 
The  format  for  this  name  is  A7  (left-justified)  and  it  may  contain 
any  series  of  seven  or  less  alphanumeric  characters,  including 
spaces.  After  reading  this  name,  the  program  generates  the  random 
variables  and  writes  them  out  to  a  disk  file.  At  this  point  the 
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program  writes  out  the  file  name  and  a  summary  of  namelist 
ables  on  unit  6  and  then  looks  to  unit  5  for  more  namelist 
Piles  are  generated  in  this  manner  until  an  end-of-file  is 
on  unit  5,  at  which  point  execution  terminates. 


vari- 
input . 
detected 


NRMGEN  —  Gaussian  Random  Number  Generator 

This  program  calls  the  IMSL  subroutine  GGNPN  to  generate 
independent,  zero-mean,  unit-variance,  Gaussian  random  variables. 
The  loader  must  be  instructed  to  search  the  IKSL  library  for 
external  references  when  NRMGEN  is  loaded. 

The  random  variables  are  written  out  to  a  disk  file  in  1000- 
variable  segments.  The  operator  can  instruct  the  program  to 
square  each  variable  before  writing  it  to  disk.  The  first  file 
is  written  on  logical  I/O  unit  20.  Succeeding  files  are  written 
on  unit  21,  unit  22,  and  so  on.  Output  summarizing  the  content 
of  each  file  is  written  to  unit  6. 

The  following  variables  comprise  the  namelist  for  NRMGEN: 

31 

DSEED  -  a  real  number  between  0  and  2  exclusive 
that  is  used  to  "seed"  the  generation 
process.  This  number  should  only  be 
specified  in  the  first  of  the  series  of 
namelist  input  for  any  particular  execu¬ 
tion.  Otherwise,  the  random  variables 
may  not  be  independent  from  file  to  file. 

This  is  because  each  call  to  GGNPM  uses  the 
seed  returned  from  the  last  call  and 
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generates  a  new  seed  to  be  used  by  the 
succeeding  call.  Continuity  of  this  chain 
insures  independence.  After  generation  of 
each  file,  the  last  value  of  DSEED  is  printed 
out  in  the  file  summary  at  unit  6.  On  future 
executions  of  NRMGEN,  DSEED  should  be  set 
equal  to  the  last  output  value  of  the  pre¬ 
vious  execution  to  insure  that  the  new  files 
are  independent  of  the  previously  generated 

ones.  The  subroutine  GGNPM  should  not  repeat 
31 

until  after  2  variables. 

NUM  —  an  integer  that  specifies  the  number  of  1000- 
variate  segments  to  be  written  to  the  disk 
file. 

SQ  —  a  logical  variable  that  determines  whether 
the  variables  are  squared  before  being 
written  to  disk.  This  variable  should  be 
set  equal  to  ".true.”  if  the  variables  are 
to  be  squared.  Otherwise  it  should  be  left 
out  of  the  namelist. 

EXPGEN  -  Exponential  Random  Number  Generator 

This  progarm  calls  the  IMSL  subroutine  GGUBS  to  generate 
independent  random  variables  uniformaly  distributed  between 
0  and  1.  These  variables  are  converted  to  independent,  unit 
mean  exponential  variables  through  the  transformation  E  *  EXP(u) 
the  loader  must  be  instructed  to  search  the  IMSL  library  for 
external  references  when  EXPGEN  is  loaded. 
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The  random  variables  are  written  out  to  a  disk  file  in  1000- 
variable  segments.  The  first  file  is  written  on  logical  I/O  unit 
20.  Succeeding  files  are  written  on  unit  21,  unit  22,  and  so  on. 

The  following  variables  comprise  the  namelist  input  for 
EXPGEN: 

•31 

DSEED  —  a  real  number  between  0  and  2  that  is  used 
to  "seed"  the  generation  process.  Experience 
has  shown  that  this  number  should  be  greater 
than  one  billion,  to  insure  independence  across 
files,  this  number  should  only  be  specified  in 
the  first  of  the  series  of  namelist  input  for 
any  particular  execution  of  EXPGEK'  (see  the 
instructions  for  DSEED  in  NRMGEN). 

NUM  -  an  integer  that  specifies  the  number  of  1000- 

variate  segments  to  be  written  to  the  disk  file. 

CHI6EN  -  Chi-Squared  Random  Number  Generator 

This  program  reads  previously-generated,  independent  Gaussian 
variables  as  input  and  generates  Chi-squared  random  variables  as 
output.  The  block  diagram  for  the  generation  process  is  shown 
in  Figure  2.4.  The  Gaussian  variables  are  input  from  a  file 
attached  to  logical  I/O  unit  30.  This  file  is  rewound  and  read 
in  again  from  the  beginning  for  every  succeeding  file  of  Cni- 
squareds  that  is  generated.  The  Chi-squareds  are  written  out  to 
a  disk  file  in  1000-variable  segments.  The  operator  has  the 
option  of  writing  the  square  root  of  every  variable  out  to  disk 
rather  than  the  variable  itself.  The  first  file  of  Chi-squareds 
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Is  written  out  to  logical  I/O  unit  20.  Succeeding  files  are 
written  out  to  unit  21,  unit  22,  and  so  on. 

The  following  variables  comprise  the  namelist  input  for 
CHI GEN : 

NDF  —  the  number  of  degrees  of  freedom  of  the 
Chi-squared  series. 

NRA  -  the  length  of  the  running  average  in  the 
generation  process. 

NUK  —  the  number  of  1000-variable  segments  to  be 
generated 

SQ  —  a  logical  variable  that  determines  whether 
the  square  root  of  each  Chi-squared  or  the 
Chi-squared  Itself  is  written  to  disk.  This 
variable  should  be  set  to  ".true."  if  square 
roots  are  desired.  Otherwise,  SQ  should  be 
left  out  of  the  namelist  input. 

B-2.  Multi-Channel  Processor  Programs 

These  models  simulate  the  passive  sonar  detection  process. 
Each  model  contains  3  channels.  Signals  representing  noise  are 
always  input  at  all  three  channels.  The  middle  channel  may 
also  contain  a  genuine  signal.  The  channel  inputs  are  processed 
to  detect  signals  present  at  the  center  channel.  The  outside 
channels  are  needed  because  the  detection  procedure  works  by 
comparing  the  center  channel  to  the  outside  ones. 
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In  the  interest  of  computational  efficiency,  the  simulations 
were  divided  into  four  categories  and  a  specialized  model  was  built 
for  each  category.  The  simulations  were  divided  according  to 
whether  they  were  used  to  estimate  detection  or  false  alarm  prob¬ 
abilities,  and  whether  the  power  envelopes  of  the  noise  series 
were  fully  or  only  partially  correlated  across  each  channel. 

Each  program  is  written  in  Fortran.  The  programs  receive 
the  simulation  instructions  and  parameters  (length  of  simulation, 
degrees  of  freedom  of  the  noise  envelope,  etc.)  through  the  name- 
list  FACT  at  logical  I/O  unit  5*  Each  program  reads  from  unit 
5  until  the  namelist  is  satisfied  (see  a  Fortran  manual  for  details 
on  how  to  satisfy  a  namelist)  at  which  point  simulation  is  begun. 
After  the  first  simulation  the  program  reads  again  from  unit  5, 
and  performs  another  simulation.  The  program  ends  when  an  end  of 
file  is  encountered  at  unit  5- 

MCHCFA  —  Common  Noise  Power  Envelope  False  Alarm  Model 

This  model  is  designed  to  estimate  the  false  alarm  probability 
of  the  multichannel  processor  when  the  power  envelope  of  the  noise 
is  common  to  all  three  channels.  The  block  diagram  for  this  model 
is  shown  in  Figure  5*1.  To  save  CPU  time,  the  squared,  zero-mean, 
unit-variance,  Gaussian  n. ise  series  at  the  two  outside  channels 
are  replaced  by  a  series  of  independent  unit-mean,  exponential 
variables  before  multiplication  by  the  power  envelope.  The  power 
envelope  is  always  a  Chi-squared  series  generated  by  the  CHI GEN . 

The  input  to  the  center  channel  is  always  a  series  of  independent, 
squared,  Gaussian  variates. 
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The  following  variables  comprise  the  namelist  input  for 
MCHCFA: 

NCHIF  -  the  logical  unit  number  of  the  file  that  contains 
the  power  envelope.  The  default  for  this  variable 
is  20. 

NDF  —  the  number  of  degrees  of  freedom  in  the  Chi- 

squared  random  variables  in  the  power  envelope. 

NPRA  —  the  length  of  the  running  average  in  the  power 
envelope . 

NRECRA  —  the  length  of  the  rectifier  running  average  in 
the  MC  processor. 

NSIM  -  the  number  (in  thousands)  of  independent  obser¬ 
vations  of  the  rectifier  output. 

NC  —  the  number  of  different  values  for  parameter  C 
in  this  simulation. 

C  —  this  is  a  vector  (maximum  length  10)  of  values 
for  the  threshold  parameter  that  is  used  in¬ 
ternally  by  the  processor.  False  alarm  estimates 
are  generated  for  each  value  of  C  during  any 
simulation. 

The  Gaussian  input  for  MCHDFA  is  read  at  unit  30.  The 
exponential  input  is  read  at  unit  31* 
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OUTPUT 


The  probability  of  false  alarms  for  each  value  of  C  consti¬ 
tutes  the  most  important  part  of  the  output.  This  probability 
is  merely  the  number  of  detections  divided  by  the  number  of 
independent  observations.  The  names  of  the  files  used  as 
input,  and  a  summary  of  the  ‘'facts"  namelist  are  also  output. 

MCtiCDT  —  Common  Noise  Power  Envelope  Detection  Model 

This  model  is  designed  to  estimate  the  detection  probability 
of  the  processor  when  the  power  envelope  of  the  noise  is  common 
to  all  channels.  The  block  diagram  for  this  model  is  shown  in 
figure  6.1.  This  model  is  exactly  the  same  as  the  previous  one 
except  that  an  additional  series  is  added  to  the  center  channel 
to  simulate  a  signal.  This  additional  signal  is  multiplied  by  a 
power  envelope  that  is  independent  of  the  noise  power  envelope. 
The  signal  power  envelope  is  generated  by  the  Chi-squared  program 
and  is  the  square  root  of  a  Chi-squared  series.  The  two  series 
at  the  bottom  left  of  Figure  6.1  represent  this  additional  input. 

A  list  of  the  variables  that  comprise  the  namelist  follows: 

NCHIS  —  The  logical  unit  number  of  the  file  that 

contains  the  power  envelope  of  the  signal. 

The  default  is  40. 

NCHIN  -  the  logical  unit  number  of  the  file  that 
contains  the  power  envelope  of  the  noise. 

The  default  of  this  variable  is  20. 
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NDFS  —  the  number  of  degrees  of  freedom  in  the  running 
average  of  the  signal  envelope. 

NDFH  —  the  running  average  in  the  noise  envelope. 

NPRAS  —  the  running  average  in  the  signal  envelope 

NRECRA  —  the  length  of  the  running  average  in  the 
rectifier . 

NSIK  —  the  number  (in  thousands)  of  independent 
observations  of  the  rectifier  output. 

NK  —  the  number  of  different  values  for  the 
S/N  ration  in  this  simulation. 

SN  -  a  vector  (maximum  length  10)  of  signal-to~ 
noise  ratios.  A  detection  probability  is 
given  for  each  value  of  S/N. 

C  —  the  threshold  parameter ...  used  internally 
by  the  processor. 


OUTPUT 


The  detection  probability  for  each  value  of  S/N  is 
output  along  with  the  names  of  the  files  used  as  input  and  a 
summary  of  the  "facts"  namelist. 
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MCHPDT  -  Partially-Correlated  Noise  Power  Envelope  Detection  Model 

This  model  is  designed  to  estimate  the  detection  probability 
of  the  processor  when  the  power  envelopes  of  the  noise  inputs  are 
partially  correlated  across  the  3  channels.  The  block  diagram 
for  this  model  is  shown  in  Figure  7.2.  Note  that  the  outside 
channels  cannot  be  combined  for  this  model  because  of  their 
different  power  envelopes. 

The  generation  procedure  is  internal  to  MCHDDT  and  not  a 
part  of  CHIGEN.  Only  one  series  of  independent,  zero-mean, 
unit-variance,  Gaussian  variates  is  required  as  input  by  this 
procedure.  As  with  MCHCDT,  the  power  envelope  for  the  signal  is 
Chi-squared  and  independent  of  the  noise  power  envelope. 

A  list  of  the  variables  that  comprise  the  namelist  input 
follows : 

NCHIS  —  the  logical  I/O  unit  number  of  the  file  that 
contains  the  signal  power  envelope  series. 

NCHIN  —  the  logical  I/O  unit  number  of  the  file  that 
contains  the  series  of  independent  Gaussian 
variates  that  drive  the  partially  correlated 
Chi-squared  generator. 

NDFS  —  the  number  of  degrees  of  freedom  of  the  power 
envelope  of  the  signal. 

NDFN  —  the  number  of  degrees  of  freedom  of  the  power 
envelope  of  the  noise. 
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NFRAS  —  the  length  of  the  running  average  in  the 
power  envelope  of  the  signal. 

NPRAN  —  the  length  of  the  running  average  in  the  power 
envelope  of  the  noise. 

NRECRA  —  the  length  of  the  running  average  following 
the  rectifier. 

NSIM  —  the  number  (in  thousands)  of  independent  obser¬ 
vations  of  rectifier  output.  The  time  steps 
between  observations  to  achieve  independence 
is  computed  internally  by  the  program  is 

NRECRA  +  ma  (NPRAS,  NPRAN). 

A 

NK  —  the  number  of  different  values  for  SN  in 
this  simulation. 

SN  —  a  vector  S/N  values  (maximum  length  10). 

C  —  the  threshold  parameter  used  internally 
by  the  processor. 


OUTPUT 

The  detection  probability  for  each  value  of  SN  is  output 
along  with  the  names  of  files  used  as  input  and  a  summary  of 
the  "facts"  namelist. 
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MCHPFA  —  Partially  Correlated  Noise  Power  Envelope  False  Alarm 
Model 


This  model  is  designed  to  estimate  the  false  alarm  probability 
of  the  processor  when  the  power  envelope  of  the  noise  is  partially 
correlated  across  the  3  channels.  The  only  difference  between 
this  model  and  MCKPDT  is  that  there  is  no  signal  series  input 
to  the  center  channel  in  this  model.  No  additional  block  diagram 
is  shown  for  this  model. 

A  list  of  the  input  variables  in  the  namelist  for  this 
program  follows: 

NCHIN  —  The  logical  I/O  unit  of  the  file  that  contains 
the  series  of  independent  Gaussians  that  drive 
the  partially  correlated  Chi-squared  genera¬ 
tor. 

NDFN  —  the  number  of  degrees  of  freedom  of  the  power 
envelope  of  the  signal. 

NPRAN  —  the  length  of  the  running  average  in  the  power 
envelope  of  the  noise. 

NRECRA  —  the  length  of  the  running  average  following 
rectification. 

NK  —  the  number  of  different  values  for  the 
parameter  C. 

C  —  a  vector  of  values  for  the  threshold  parameter 
(maximum  length  10). 


i 
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OUTPUT 

The  false  alarm  probability  for  each  value  of  C  is  output 
along  with  the  names  of  the  files  used  as  input  and  a  summary 
of  the  variables  in  the  "facts"  namelist. 
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